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Inhomogeneous Tsallis distributions in the HMF model 
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Abstract. We study the maximization of the Tsallis functional at fixed mass and energy in the HMF 
model. We give a thermodynamical and a dynamical interpretation of this variational principle. This leads 
to q-distributions known as stellar polytropes in astrophysics. We study phase transitions between spatially 
homogeneous and spatially inhomogeneous equilibrium states. We show that there exists a particular index 
q c = 3 playing the role of a canonical tricritical point separating first and second order phase transitions in 
the canonical ensemble and marking the occurence of a negative specific heat region in the microcanonical 
ensemble. We apply our results to the situation considered by Antoni & Ruffo [Phys. Rev. E 52, 2361 (1995)] 
and show that the anomaly displayed on their caloric curve can be explained naturally by assuming that, 
in this region, the QSSs are polytropes with critical index q c = 3. We qualitatively justify the occurrence 
of polytropic (Tsallis) distributions with compact support in terms of incomplete relaxation and inefficient 
mixing (non-ergodicity) . Our paper provides an exhaustive study of polytropic distributions in the HMF 
model and the first plausible explanation of the surprising result observed numerically by Antoni & Ruffo 
(1995). In the course of our analysis, we also report an interesting situation where the caloric curve presents 
both microcanonical first and second order phase transitions. 

PACS. 5.20.-y Classical statistical mechanics - 05. 45. -a Nonlinear dynamics and chaos - 05.20.Dd Kinetic 
theory - 64.60.De Statistical mechanics of model systems 



1 Introduction 

Systems with long-range interactions are numerous in 
nature. Some examples include self-gravitating systems 
(galaxies), two-dimensional turbulence (vortices), chemo- 
taxis of bacterial populations (clusters) and some mod- 
els in plasma physics [Tj. These systems are fascinating 
because they present striking features that are absent 
in systems with short-range interactions such as nega- 
tive specific heats in the microcanonical ensemble, nu- 
merous types of phase transitions, ensembles inequiv- 
alence, unusual thermodynamic limit, violent collision- 
less relaxation, long-lived quasi stationary states (QSS), 
non-Boltzmannian distributions, out-of-equilibrium phase 
transitions, re-entrant phases, non-ergodic behavior, slow 
collisional relaxation, dynamical phase transitions, alge- 
braic decay of the correlation functions... We refer to |2l3j 
for some recent reviews on the subject. 

In order to understand these strange properties in a 
simple setting, a toy model of systems with long-range 
interactions has been actively studied. It consists of N 
particles moving on a ring and interacting via a cosine po- 
tential. This model has been introduced by many authors 
[4-8] at about the same period with different motivations 
(see a short history in [5])- Let us mention for example that 
Pichon [7] introduced this model to explain the formation 



of bars in disk galaxies, giving it a physical application. 
This model is now known as the Hamiltonian Mean Field 
(HMF) model. This abbreviation was coined by Antoni & 
Ruffo [5] and it stood at that time either for Hamiltonian 
or Heisenberg Mean Field model since the HMF model 
can be viewed as a mean field XY model with long-range 
interactions (i.e. not restricted to the nearest neighbors). 
In fact, this model was first introduced much earlier (this 
is not well-known) by Messer & Spohn [TU] and called the 
cosine modeQ- They rigorously studied phase transitions 
in the canonical ensemble by evaluating the free energy 
and exhibited a second order phase transition between a 
homogeneous phase and a clustered phase below a critical 
temperature T c = 1/2. 

In their seminal paper, Antoni & Ruffo [5] studied the 
statistical mechanics of this model in the canonical en- 
semble (CE) directly from the partition function and per- 
formed A^-body simulations in the microcanonical ensem- 
ble (MCE). They started from a waterbagQ initial condi- 



was sug- 



1 These authors mention that this "cosine model" 
gested by G. Battle [IT] . 

2 A waterbag distribution corresponds to a uniform distribu- 
tion function f(8, v) = fo in the domain [— 9 m , 9 m ] x [—v m ,v m ] 
surrounded by "vaccum" f(9, v) = 0. It can have different mag- 
netization < Mq < 1. 
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tion with magnetization Mq = M(0) = 1 and plotted in 
their Fig. 4 the caloric curve giving the averaged kinetic 
temperature T = 2Eki n /N as a function of the energy U . 
They compared their numerical curve with the Boltzmann 
prediction of statistical equilibrium in the canonical en- 
semble. They found a good agreement at high and low en- 
ergies. However, close to the critical energy U c = 3/4 (cor- 
responding to T c = 1/2), the results differ from the canon- 
ical prediction. In particular, the system bifurcates from 
the homogeneous branch at a smaller energy U ~ 0.6 — 0.7 
and the caloric curve T(U) displays a region of negative 
specific heats. Antoni & Ruffo [Hj interpreted this result 
cither as (i) a nonequilibrium effect, or (ii) as a manifesta- 
tion of ensembles inequivalence due to the non-additivity 
of the energy. However, it is clear from the previous work 
of Inagaki [5] and the application of the Poincare theorem 
that the ensembles are equivalent for the HMF model (as 
confirmed later by various methods [121912] ). Therefore, 
this negative specific heats region is not a manifestation 
of ensembles inequivalence but rather a nonequilibrium 
effect. 

Latora et al. |13ll4j again performed microcanoni- 
cal simulations of the HMF model. They confirmed the 
"anomaly" (negative specific heats region) reported by 
Antoni & Ruffo [8] and observed many other anomalies in 
that region such as non-Gaussian distributions, anomalous 
diffusion, Levy walks and dynamical correlations in phase- 
space. Furthermore, they showed that the thermodynamic 
limit — > +oo and the infinite time limit t — > +oo do 
not commute. They evidenced two regimes in the dynam- 
ics. On a short timescale, of the order of the dynamical 
time irj ~ 1, the system reaches a quasistationary state 
(QSS). Then, on a much longer timescale t re i ax (N), the 
system relaxes towards the Boltzmann distribution of sta- 
tistical equilibrium. They showed that the relaxation time 
increases rapidly (algebraically) with the number of parti- 
cles N so that, at the thermodynamic limit N —> +oo, the 
system remains permanently in the QSS. It is thus clear 
from this study that the different anomalies mentioned 
above characterize the QSS, not the statistical equilib- 
rium state that is reached much later. In particular, the 
results of Antoni & Ruffo [8] and Latora et al. [14] show 
that the QSS is non-Boltzmannian in the region close to 
the critical energy U c . Latora et al. [T3] thus proposed to 
describe the QSS in terms of Tsallis [15] generalized ther- 
modynamics leading to q-distributions. Note that there is 
no reason why the QSS should be Boltzmannian (in the 
usual sense) since it is an out-of-equilibrium structure. 
Therefore, the comparison of the numerical caloric curve 
with the Boltzmann equilibrium caloric curve T(U) is not 
justified a priori. 

After the conference in Les Houches in 2002, and in- 
spired by the results in astrophysics and 2D turbulence 
presented by one of the authors [T5], several groups of 
researchers [17I9I18I1T)] started to interpret the QSSs ob- 
served in the HMF model in terms of Lynden-Bell's [5D] 
statistical theory of violent relaxation based on the Vlasov 
equation. This is a fully predictive theory based on stan- 
dard thermodynamics but taking into account the speci- 



ficities of the Vlasov equation (Casimir constraints). Ap- 
plying the Lynden-Bell theory to the HMF model (for 
waterbag initial conditions), an out-of-cquilibrium phase 
transition was discovered in [18I19J . For a given value 
of the energy U, there exists a critical value of the 
initial magnetization {Mo) C rit(U) such that: for Mq < 
(Mo) cr it(U) the stable Lynden-Bell distribution (i.e. most 
probable state) is homogeneous (non-magnetized) and 
for Mo > (Mo)erit(U) the stable Lynden-Bell distribu- 
tion is inhomogeneous (magnetized). For U = 0.69, the 
critical magnetization is (Mo) cr i t = 0.897 [18119] . More 
generally, there exists a critical line {M ) crit (U) in the 
phase diagram separating homogeneous and inhomoge- 
neous Lynden-Bell distributions. It was found later [21122] 
that the system displays first and second order phase tran- 
sitions separated by a tricritical point. There is also an 
interesting phenomenon of phase reentrance in the (/o, U) 
plane predicted in [18] and numerically confirmed in |23) . 
Coming back to the specific value U = 0.69, direct nu- 
merical simulations of the HMF model for Mq < (Mo) cr a 
have shown that the Lynden-Bell prediction works fairly 
well [T5]. This agreement is remarkable since there is no 
fitting parameter in the theory. This led the authors of 
[T9] to argue that "Lynden-Bell's theory explains quasista- 
tionary states in the HMF model" . A controversy started 
when some of these authors [1712412] concluded that "the 
approach of Tsallis is unsuccessful to explain QSSs" . 

However, caution was made by one author [25118126] 
who argued that the Lynden-Bell theory does not explain 
everything. Indeed, for initial magnetization Mq = 1, the 
system is in the non-degenerate limit so that the Lynden- 
Bell entropy reduces to the Boltzmann entropy (with a 
different interpretation). Therefore, in this limit case, the 
Lynden-Bell theory leads exactly to the same prediction 
as the usual Boltzmann statistical theory but, of course, 
for a completely different reason. This observation led to 
a re-interpretation [55] of the caloric curve T(U) obtained 
by Antoni & Ruffo [8] and Latora al. [14]. In this curve, 
the theoretical line should not be interpreted as the Boltz- 
mann statistical equilibrium state but as the Lynden-Bell 
statistical equilibrium stat^E They turn out to coincide in 
the case Mq = 1 but this is essentially coincidental. With 
this new interpretation [26) . the comparison reveals that 
the Lynden-Bell theory works well for large and low ener- 
gies but that it fails for energies close to the critical energy. 
Therefore, the Lynden-Bell theory does not explain the ob- 
servations in the range [0.5, U c ]. Chavanis [26] interpreted 
this disagreement as a result of incomplete relaxation. In- 
deed, it was emphasized by Lynden-Bell [20] himself that 
his approach implicitly assumes that the system "mixes 
efficiently" so that the ergodicity hypothesis which sus- 
tains his statistical theory applies. However, it has been 
observed in many cases of violent relaxation that the sys- 
tem does not mix efficiently so that the Lynden-Bell pre- 



3 As mentioned above, there is no reason to compare the 
QSS with the Boltzmann prediction that corresponds to the 
collisional regime reached for t —¥ +oo. By contrast, it is fully 
relevant to compare the observed QSS with the Lynden-Bell 
prediction that applies to the collisionless regime. 
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diction fails (see various examples quoted in [25126] L The 
qualitative reason is easy to understand |26| . Since vio- 
lent relaxation is a purely inertial process (no collision), 
mixing is due to the fluctuations of the mean field poten- 
tial caused by the fluctuations of the distribution function 
itself [2DJ. However, as the system approaches metaequilib- 
rium (QSS), the fluctuations of the distribution function 
are less and less efficient (by definition!) and the system 
can be trapped in a QSS (a steady solution of the Vlasov 
equation on the coarse-grained scale) that is not the most 
mixed state. This is what happens in astrophysics (ellip- 
tical galaxies are not described by Lynden-Bell's distribu- 
tion that has infinite mass [27] ) and in certain situations 
of 2D turbulence [28-30]. The results of Antoni & Ruffo [S] 
and Latora et al. [14j reveal that the same phenomenon oc- 
curs for the HMF model close to the critical energ jfl This 
is also visible on Fig. 8 of Bachelard et al. [3T] . There is a 
huge region in the top right of the phase diagram where 
the Lyndcn-Bcll theory does not work. This concerns in 
particular the point U = 0.69 and Mq = 1, as anticipated 
in [18126] . It is precisely this "no-man's land" region that 
Tsallis and coworkers have investigated [3] • In this region, 
standard statistical mechanics (i.e. Lynden-Bell's theory) 
does not seem to directly appl}0. 

In their early work, Latora et al. [2] tried to fit the 
QSS by a q-distribution. They considered a distribution 
with q = — 5 < 1 (in our notations) leading to power- 
law tails and introduced by hand an additional cut-off 
to make the distribution normalizable. This procedure is 
very ad hoc. Furthermore, even if we accept it, we can ar- 
gue that it does not provide an impressive fit of the QSS. 
Recently, Campa et al. [34j have performed new simula- 
tions for initial magnetizations Mo = and Mo = 1 and 
found that the QSS is very well-fitted by a semi- ellipscB 
Similar results were obtained earlier by Yamaguchi et al. 
[T7] for the Mq — case. They claimed that the QSS is 

4 Latora et al. [14] find that the Largest Lyapunov exponent 
for the QSS tends to zero. In this sense, mixing is negligible and 
one expects anomalies in the relaxation process. This is fully 
consistent with the idea of incomplete relaxation and inefficient 
mixing introduced by Lynden-Bell [20] and further discussed 
by Chavanis [25I1&I26] . 

5 Quoting Einstein [32] and Cohen [33] . Latora et al. [14] 
argued that standard statistical mechanics fails when the dy- 
namics plays a nontrivial role (e.g. long-range correlations or 
fractal structures in phase space). In our point of view, this 
is a correct interpretation although standard thermodynamics 
should refer here to Lynden-Bell's theory, which is the proper 
Boltzmann approach applied to the Vlasov equation. This sub- 
tlety is not addressed in the paper of Latora et al. [14] since 
they did not know the theory of Lynden-Bell at that time. Yet, 
their general comment can be applied to Lynden-Bell's theory 
as well: if the dynamics is nontrivial and the system does not 
mix well, standard (Lynden-Bell) statistical mechanics fails. 

6 For Mo = 1, this differs from the numerical results of La- 
tora et al. Q3] . However, Campa et al. [M] show that the ordi- 
nary waterbag initial condition leads to the presence of large 
sample to sample fluctuations so that many averages are nec- 
essary. They proposed to use isotropic waterbag distributions 
to reduce the fluctuations. 



not a ^-distribution since it does not have power-law tails. 
However, Chavanis [26j remarked that a semi-ellipse is a 
^-distribution with q = 3! Since q > 1, this distribution 
has a compact support. Therefore, the numerical results 
of Campa et al. [33] and Yamaguchi et al. [T7] show that 
the system tends to select a particular Tsallis distribu- 
tion as a QSS3- Furthermore, the index q = 3 seems to 
play a particular role since Campa et al. [34] obtained the 
same index q = 3 in different situations (see their Fig. 4). 
However, until now, the reason for this particular value 
remains unknown. The fact that the QSS has a compact 
support is relatively natural in the phenomenology of in- 
complete violent relaxation. Indeed, when mixing is not 
very efficient, we expect that the high energy states are 
not sampled by the system. This leads to a confinement of 
the distribution which is a virtue of the Tsallis distribu- 
tions with q > 1 . A similar confinement was observed in a 
plasma experiment [28] and a good fit was obtained with 
a g-distribution with index q = 2 (in our notations) |29) . 
This confinement was justified by a lack of ergodicity in 
the system [30] . One can therefore interpret the Tsallis dis- 
tributions as an attempt to take into account incomplete 
relaxation and non-ergodicity in systems with long-range 
interactions. In this interpretation, the index q could be a 
measure of mixing. If we assume that the system mixes 
efficiently, then q = 1 and we get the Lynden-Bell theory. 
If the system does not mix well, the Lynden-Bell theory 
fails and q ^ 1. Since the value of the q parameter de- 
pends on the efficiency of mixing (which is not known a 
priori), it appears difficult to determine its value from first 
principles. Furthermore, its value can change from case to 
case since the degree of mixing can vary depending on the 
initial condition (some systems can mix well and others 
less). Finally, we can argue that the Tsallis entropy is just 
one generalized entropy among many others and that it 
may not be universal [36] . It may just describe a special 
type of non-ergodic behavior but not all of them. In fact, 
non-ergodic effects can be so complicated that it is hard 
to believe that they can be encapsulated in a simple func- 
tional such as the Tsallis functional or any other |3QI25j . 
Nevertheless, we must recognize that some QSSs can be 
well-fitted by ^-distributions. Even more strikingly, follow- 
ing the suggestion of [26] , Campa et al. [37] demonstrated 
numerically that, during the collisional regime, the time 
dependent distribution f(v,t) is still very well- fitted by 
g(i)-distributions with a time-dependent index. When the 
index reaches a critical value q C rit{U) (predicted by the 
theory [26137] ). the distribution function becomes Vlasov 



7 Let us be more precise: (i) For Mo = 1, an isotropic wa- 
terbag distribution violently relaxes towards a q-distribution 
with q = 3 [33]; (ii) For Mo = 0, a waterbag distribution is 
Vlasov stable and does not undergo violent relaxation (it is al- 
ready the Lynden-Bell state). However, in the collisional regime 
(due to finite N effects), it becomes a (/-distribution with q — 3 
[17134] : (iii) for intermediate values of Mo, the system violently 
relaxes towards the Lynden-Bell distribution [19 34 . 

8 This interpretation was proposed by one of the authors in 
several papers 35 18 26 and it may be more accurate than the 
usual interpretation: u q is a measure of nonextensivity" [3]. 
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unstable and a dynamical phase transition from the ho- 
mogeneous phase (non-magnetized) to the inhomogeneous 
phase (magnetized) is triggered. This explains previous 
observations on the evolution of the magnetization M(t) 
[14I17I34] . Interestingly, a very similar behavior has been 
found by Taruya & Sakagami [35] for self-gravitating sys- 
tem^- 

In view of these results, it is interesting to study 
in more detail the structure and the stability of q- 
distributions. Note that g-distributions correspond to 
what have been called stellar polytropes in astrophysics 
[27] . They were introduced long ago by Eddington [3D] 
as particular stationary solutions of the Vlasov equation. 
They were used to construct simple self-consistent mathe- 
matical models of galaxies. At some time, they were found 
to provide a reasonable fit of some observed star clusters, 
the so-called Plummer [3T| model. Improved observation 
of globular clusters and galaxies showed that the fit is not 
perfect and more realistic models have been introduced 
since then [57]. However, stellar polytropes are still im- 
portant for historical reasons and for their mathematical 
simplicity. The stability of polytropic distributions is an 
old problem in stellar dynamics [27J. It has been recon- 
sidered recently, for box-confined systems, by Taruya & 
Sakagami [42-44] in the framework of Tsallis generalized 
thermodynamics and by Chavanis et al. [45 35 46 147139] in 
relation to their nonlinear dynamical stability with respect 
to the Euler and Vlasov equations. It is therefore interest- 
ing to extend these studies to the case of the HMF model. 
For the moment, only spatially homogeneous polytropic 
distributions have been considered |9l48j . In the present 
paper, we extend the analysis to spatially inhomogeneous 
polytropes. Specifically, we study the maximization of the 
Tsallis functional at fixed mass and energy and plot the 
corresponding caloric curves. We give a thermodynamical 
and a dynamical interpretation of this variational princi- 
ple. We find the existence of a critical polytropic index 
q c = 3 (where results are analytical) which plays the role 
of a canonical tricritical point separating first and sec- 
ond order phase transitions in the canonical ensemble and 
marking the onset of negative specific heats in the mi- 
crocanonical ensemble (there also exists a microcanonical 
tricritical point at q ~ 16.9 and a microcanonical critical 
point at q ~ 6.55). Interestingly, this critical value q c = 3 
turns out to be the one observed by Campa al. [33] in their 
numerical simulations. Then, we apply our theory to the 
situation considered by Antoni & Ruffo [SJ. We find that 
the structure of their numerical caloric curve T{U) close 
to the critical energy can be explained naturally if we as- 
sume that the QSSs in this region are polytropes with the 
critical index q c = 3 (this result is specifically discussed in 
Sec. 17. 2p . This yields a transition energy U' c = 5/8 lower 
than U c = 3/4 and a region of negative specific heats in 
the curve T(U) that are in qualitative agreement with the 
numerical results (the agreement could be improved by al- 



There is, however, a crucial difference in the interpretation. 
Taruya & Sakagami [35J interpret the instability in terms of 
Tsallis generalized thermodynamics while Campa et al. |37| 
interpret it in terms of Vlasov dynamical instability [39] . 



lowing q to deviate slightly from the critical value q c = 3). 
These behaviors had never been explained since the orig- 
inal paper of Antoni & Ruffo [5J- We provide a plausible 
explanation in terms of Tsallis (polytropic) distributions. 
We qualitatively justify the occurence of polytropic dis- 
tributions with a compact support in terms of incomplete 
relaxation. Furthermore, the index selected by the sys- 
tem appears to be the critical one (or close to it). This 
is the first time that a sort of prediction of the q index is 
made in that context. However, our approach does not ex- 
plain everything and just opens a direction of research. In- 
deed, one has first to confirm the results by more detailed 
comparisons with numerical simulations and, in case of 
agreement, try to understand why critical polytropes are 
selected by the system and if this is a general feature. 

2 Interpretations of the Tsallis functionals 

In order to motivate our study of polytropes, we shall first 
recall different interpretations of the Tsallis functionals 
that have been given previously by one of the authors 
[47149] . 

2.1 Thermodynamical interpretation 

There is no reason to describe the QSSs that emerge in 
Hamiltonian systems with long-range interactions in terms 
of Boltzmann statistical mechanics since they correspond 
to out-of-equilibrium structures (in the usual sense). The 
QSSs are formed during the collisionless regime while 
Boltzmann statistical equilibrium is reached on a much 
longer timescale at the end of the collisional regime. 
Fundamentally, the QSSs are stable steady states of the 
Vlasov equation (on a coarse-grained scale) and they 
should be described by Lynden-Bell's statistical mechan- 
ics. However, Lynden-Bell's theory, as any statistical the- 
ory, assumes efficient mixing and ergodicity. If the system 
does not mix well, the QSS will differ from the Lynden- 
Bcll prediction (by definition). If we want to apply Tsallis 
generalized thermodynamics to that context, in order to 
take into account non ergodic effects and incomplete mix- 
ing, we must modify the Lynden-Bcll entropy which is the 
Boltzmann entropy associated to the Vlasov equation^- 
Therefore, the proper g-entropy to consider is 

S q [p] = — ^ J (p 9 - p) dndrdv. (1) 

This entropy applies to the distribution p(r, v, ij) of phase 
levels T] (see [35] for details) so as to take into account the 

10 Latora et al. [14] did not know the Lynden-Bell theory 
and proposed to replace the usual Boltzmann entropy S[f] = 
- / /ln/drdv by the g-entropy S q [f] = / (/« -/) drdv. 

However, as argued long ago in [3D], this approach is in general 
incorrect since it does not take into account the constraints 
of the Vlasov equation. Furthermore, since S q [f] is rigorously 
conserved by the Vlasov equation (it is a particular Casimir), 
there is no thermodynamical reason to maximize it. 
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constraints of the Vlasov equation (Casimirs). For q = 1 
(efficient mixing), we recover the Lynden-Bell entropy. For 
q 5^ 1, this functional could describe incomplete violent re- 
laxation and non-crgodic effects in the framework of Tsal- 
lis thermodynamics. Now, for two-levels initial conditions 
and in the non degenerate limit (this corresponds to the 
waterbag model with initial magnetization Mq = 1 in the 
HMF model), the Lynden-Bell entropy takes a form sim- 
ilar to the Boltzmann entropy and the constraints reduce 
to the mass and the energy. In that case, the generalized 
Lynden-Bell entropy (HJ) reduces to 



S q [f] 



1 



9-1 



7V 7 



no 



'/<> 



drdv. 



(2) 



This is similar to the Tsallis entropy considered by Latora 
et al. [2] but it arises here for a different reason: note 
in particular the bar on / (coarse-grained distribution) 
and the presence of the initial phase level 770- Note that 
S q [f] is expected to increase while S q [f] is conserved. In 
this thermodynamical approach, the QSS is obtained by 
maximizing the Tsallis entropy at fixed mass and energy 
(microcanonical ensemble). This is expected to select the 
most probable state under some dynamical constraints (re- 
sponsible for incomplete mixing) that are implicitly taken 
into account in the form of the entropy. The canonical en- 
semble has no justification in the present context but, as 
usual, it can be useful to provide a sufficient condition of 
microcanonical stability |50j . 

Since the Tsallis entropy (Q]) includes the Lynden-Bell 
entropy as a special case (for q = 1), it can give at least 
as good, or better, results. However, we would like to em- 
phasize several limitations of Tsallis generalized thermo- 
dynamics: (i) at present, there is no theory predicting the 
value of q. This prediction is of course very difficult since 
q is related to non-ergodic effects. For the moment, q ap- 
pears as a fitting parameter measuring the degree of mix- 
ing of a system; (ii) it is not clear that all non-ergodic pro- 
cesses can be described by the Tsallis entropy. Therefore, 
this functional is probably not universal. It may, however, 
describe a certain class of non-ergodic behavioral; (iii) it 
is not clear that non-ergodic effects can be described by 
a simple entropic functional such as the Tsallis functional 
or any other. Maybe, a better approach is to consider the 
dynamics of mixing and develop kinetic theories and re- 
laxation equations of the process of violent relaxation as 
suggested in [51I25J26I39] . This kinetic approach may be 
closer to the original ideas of Einstein [35] an d Cohen [33J . 



11 It has been shown that the Tsallis entropy satisfies a lot 
of axioms satisfied by the Boltzmann entropy. This makes this 
entropy very "natural" to generalize the Boltzmann entropy. 
However, the axiomatic approach may not be the best justifi- 
cation of an entropy. The original combinatorial approach of 
Boltzmann seems to be more relevant. The Boltzmann and 
the Lynden-Bell entropies can be derived from a combinatorial 
analysis by assuming that all the microstates are equiproba- 
ble. It would be interesting to derive the Tsallis entropy from 
a combinatorial analysis by putting some constraints on the 
availability of the microstates. 



2.2 Dynamical interpretation 

We would like to give another interpretation of Tsallis 
functional that is not related to thermodynamics and that 
docs not suffer the limitations exposed above. It will show 
that the Tsallis formalism can be useful in a dynamical 
context, using a thermodynamical analogy [47139149] . 

It has been known for a long time that any distribution 
function of the form / = /(e), where e = v 2 /2 + & is the 
individual energy, is a steady state of the Vlasov equation 
[52j . In particular, ^-distributions have been introduced 
long ago by Eddington [40| in astrophysics where they 
are called stellar polytropetQ. It is also known that these 
distributions are critical points of a certain functional of 
the form 



S[f] 



C(f) drdv, 



(3) 



where C is a convex function, at fixed mass and energy. 
Furthermore, if they are maxima of this functional, they 
are nonlincarly dynamically stable with respect to the 
Vlasov equation. As far as we know, this dynamical sta- 
bility criterion was first stated by Ipser & Horwitz [53] in 
astrophysics. Furthermore, Ipser [55] introduced the func- 
tional S — — j / 1 + 1 /(" -3 / 2 ) drdv (in his notations) asso- 
ciated to stellar polytropes. This is nothing but the Tsallis 
functional. 

Let us be more precise and more general. The maxi- 
mization problem 



mzx{S[f}\E[f]=E, M[f]=M}, 



(4) 



determines a steady state of the Vlasov equation of the 
form / = /(e) with /'(e) < that is nonlinearly dynam- 
ically stable. This has been stated by Ellis et al. [55] in 
2D turbulence (for the 2D Eulcr equation) and by Ipser & 
Horwitz [51] , Tremaine et al. [57] and Chavanis [3j|] in stel- 
lar dynamics. The maximization problem ([4]) is similar to 
a condition of microcanonical stability in thermodynam- 
ics. Therefore, we can develop a thermodynamical analogy 
[49j to investigate the nonlinear dynamical stability prob- 
lem. In this analogy, S is called a pseudo entropy. Thus, 
the Tsallis functional is a particular pseudo entropy whose 
maximization at fixed mass and energy determines dis- 
tributions (polytropes) that are nonlinearly dynamically 
stable with respect to the Vlasov equation [47139] . 



The connection between g-distributions and stellar poly- 
tropes was first mentioned by Plastino & Plastino [53] . A more 
detailed discussion has been given recently by Chavanis & Sire 

13 Indeed, the second term in the Tsallis functional is a 
constant (proportional to mass) and the first term coincides 
with the Ipser functional with the notation (|18l) . As noted in 
[47] . the Tsallis functional is more convenient to take the 
limit g — > 1 (i.e. n — > +00) since it reduces to the Boltz- 
mann functional while the Ipser functional takes a trivial form. 
Therefore, the Tsallis functional allows to make a continuous 
link between isothermal (n = 00) and polytropic (n finite) dis- 
tributions using L'Hopital's rule. 
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The minimization problem 

min {F[f} = E[f] - TS[f] \ M[f] = M} . (5) 

also determines a steady state of the Vlasov equation that 
is nonlinear ly dynamically stable. In fact, (Q| and ([5]) have 
the same critical points. However, the criterion ([5]) is less 
refined than (U). Indeed, the minimization problem ([3]) is 
similar to a condition of canonical stability in thermody- 
namics. In this analogy, F is called a pscudo free energy. 
Now, it is a general result |50j that canonical stability im- 
plies microcanonical stability, but the reciprocal is wrong 
in case of ensembles inequivalence. Schematically: ([5]) => 
Q. In the present dynamical context, this means that 
steady states of the Vlasov equation that are stable ac- 
cording to the criterion (|5|) are necessarily stable accord- 
ing to the more constrained criterion (QJ. However, there 
may exist steady states of the Vlasov equation that are 
stable according to dU while they fail to satisfy ([5]). As 
shown in [35], this is the case for stellar polytropes with 
indices 3 < n < 5 in astrophysics. 

It can also be shown [35] that ([5]) is equivalent to 

Tnm{F[p]\M[p\ = M}, (6) 
p 

where 

F=±Jp$dr + JpJ P ^ dp'dv. (7) 

More precisely, a distribution function /(r, v) is solution of 
([5]) iff the corresponding density profile p(r) is solution of 
([5]), where p{p) is the equation of state determined by C(f) 
(see [49] for more details). Finally, ([6]) is clearly equivalent 
to 

mm{W\p,u]\M\p] = M}, (8) 

p,u 

where 

W = J r p J" ^^dp'dr+^ J p<Pdr + J p^-dr. (9) 

Now, it can be shown that this minimization problem de- 
termines a steady solution of the barotropic Euler equa- 
tion that is formally nonlincarly dynamically stable |27j . 
From the implication @ flU, we con- 

clude that a distribution /(r,v) is stable with respect 
to the Vlasov equation [according to Q] if the corre- 
sponding density profile p(r) is stable with respect to the 
Euler equation [according to §8$}. As shown in [39], this 
provides a nonlinear generalization of the Antonov first 
law in astrophysics: "a spherical galaxy / = /(e) with 
/'(e) < is nonlinearly dynamically stable with respect to 
the Vlasov- Poisson system if the corresponding barotropic 
star is nonlinearly dynamically stable with respect to the 
Eulcr-Poisson system" . Interestingly, this also provides a 
new interpretation [39) of this law in terms of ensembles 
inequivalence. Indeed, the Antonov first law has the same 
status as the fact that "canonical stability implies micro- 
canonical stability" in thermodynamics. 



Remark: the "microcanonical" criterion Q provides 
itself just a sufficient condition of nonlinear dynamical 
stability. There exists an even more refined criterion of 
nonlinear dynamical stability taking into account the con- 
servation of all the Casimirs (see discussion in [4"5]~). 



2.3 Generalized H-functions and selective decay 

A generalized ii-function is a functional of the coarse- 
grained distribution function /(r, v, f) of the form 

H[J] = - (c(f)dvdv, (10) 

where C is any convex function. Tremaine et al. [57] have 
shown that the generalized i7-functions increase during 
violent relaxation in the sense that H{t) > H(0) for any 
time t > where it is assumed that the initial distri- 
bution function is not mixed so that /(r, v,t = 0) = 
f(r,v,t = 00 By contrast, the energy E[J] and M[f] 
and the mass calculated with the coarse-grained distribu- 
tion function are approximately conserved. This suggests a 
phenomenological generalized selective decay principle (for 
—H): "due to mixing, the system may tend to a QSS that 
maximizes a certain _ff-function (non universal) at fixed 
mass and energy" [49p"^l In that context, the Tsallis func- 
tional H[f] = —-^i /(/ — /) drdv can be interpreted as 
particular generalized H-functions |49j . 



3 Polytropic distributions: general theory 

For any of the reasons exposed previously, we think that 
it is useful to study polytropic (Tsallis) distribution func- 
tions and investigate their stability through the variational 
problems ^ and (O . We first develop a general theory of 
polytropes following the lines of [46147] . It will be applied 
specifically to the HMF model in Sec. |3J 



3.1 Polytropic distributions in phase space 

Let us consider the Tsallis functional 

S =-^J /(/ 9 -/o 9 " 1 /)drdv. (If) 

We have introduced a constant /o in order to make the 
expression homogeneous. This constant will play no role 
in the following since the last term is proportional to the 

14 Note that the time evolution of the generalized _ff-functions 
is not necessarily monotonic (nothing is implied concerning the 
relative values of H(t) and H(t') for t,t' > 0). 

15 The close relationship between this phenomenological prin- 
ciple and the criterion of nonlinear dynamical stability @, as 
well as the limitations of this phenomenological principle, are 
discussed in detail in |25I49| . 
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mass that is conserved. Therefore, we could equally work 
with the Ipser functional 



1 



f q dvdv, 



(12) 



but the first expression allows us to make the connection 
with isothermal distributions when q — > 1. Indeed, for 
q — > 1, we recover the Boltzmann functional 



S 



/ In ( — ] drdv. 

Jo, 



(13) 



We shall consider the maximization of the Tsallis func- 
tional at fixed energy and mass 



' / y drdv + i / 



M = pdr 



(14) 



(15) 



Some interpretations of this variational problem have been 
given in Sec. [2] Here, S is either a generalized entropy 
(thcrmodynamical interpretation) or a pscudo entropy 
(dynamical interpretation). By an abuse of language, and 
to simplify the terminology, we shall call it simply an en- 
tropy. We will work in a space of dimension d since our 
formalism can have application for different systems. 

The critical points of entropy at fixed energy and mass 
are determined by the condition 



SS - (36 E - aSM = 0, 



(16) 



where (3 = 1 /T and a are Lagrange multipliers (T is the 
inverse temperature and a the chemical potential). This 
yields the ^-distributions (or polytropic distributions): 



/(r,v) 



where p = [/ 9 1 



1/(9-1) 



, (17) 

x if x > 



i g — (q—l) a] jq. The notation [x] + 
and [x]+ = if x < 0. As is customary in astrophysics, 
we define the polytropic index n by the relation 16 ! 



1 



n= 2 



9-1 



(18) 



The distribution function /(r,v) depends only on the in- 
dividual energy e = v 2 /2+<P(r), i.e. / = /(e). Therefore, it 
is a steady state of the Vlasov equation. We shall consider 
q > so that C(f) is convex and /? > so that /'(e) < 0, 
which corresponds to the physical situation. For n = d/2 
(q — > +oo), we recover the Fermi distribution [18j . For 
n — > +oo (q — > 1), we recover the isothermal distribution. 

We need to distinguish two cases depending on the sign 
of q—l. (i) For q > 1 (n > d/2), the distribution function 
can be written 



/ = A(e m - e)!" 1 , 



(19) 



16 This relation is sometimes presented as a "fundamental" 
relation relating the (/-parameter to the polytropic index n [3]- 
In our sense, this is just a definition of notations, no more. 



where we have set A = [j3(q — l)/q]'^ rr and e rn = 
qfi/{(3(q — 1)]. Such distributions have a compact sup- 
port in phase space since they vanish at e = e m . At 
a given position, the distribution function vanishes for 
v > v m (r) = v/2(e m -$(r)). (ii) For < q < 1, the 
distribution function can be written 



/ = ,4(e + e) — . 



(20) 



where we have set A = [/3(1 — q)/q]~ and e = qfi/[fi(l- 
q)}. Such distributions are defined for all velocities. At a 
given position, the distribution function behaves, for large 



velocities, as / 



- 2 /(i-9) ^ w -(d-2n), Wc shall on i y 



consider distribution functions for which the density and 
the pressure 



P= If dv, 



P=~ I fv 2 < 



are finite. This implies d/(d + 2) < q < 1 (n < — 1). 



(21) 



3.2 Polytropic equation of state 

For any distribution function of the form / = /(e), the 
density and the pressure are functions of <?(r) such that 
p = p[<P(r)] and p ~ p[<P(r)]. Eliminating <P(r) between 
these expressions, we obtain a barotropic equation of state 
p(r) = p[p(r)]. Furthermore, it is easy to see that / = 
/(e) implies the condition of hydrostatic equilibrium (see 
Appendix El : 



Vp + pW<P = 0. 



(22) 



Let us determine the equation of state corresponding to 
the polytropic distribution (jTT)) . For n > d/2, the density 
and the pressure can be expressed as (see Appendix [B]) : 

r = A S ^-*r*<*- i rW2)r r j 1 :f + n) , (23) 

r(n + 1) 



P =— (e m -<2>) 



1+1^/2-1 ^(^/2)^(1- d/2 + n) 
r(n + l) 



(24) 



where T(x) denotes the Gamma function and Sd the sur- 
face of a unit sphere in d dimensions. For n < — 1, the den- 
sity and the pressure can be expressed as (see Appendix 

p = AS d (e + ^ d/2 - l Fi r ~^^\ (25) 



P = 



^ ico+0r+ ^Et^m. (26) 



Eliminating the field <P(r) between the expressions ([23 
(jMj) and (^D-(Iini), one finds that 



p = Kp 1 , 7=1 



(27) 
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This is the classical polytropic equation of state. This is 
the reason why the distributions (|T7| are called polytropic 
distributions. Furthermore, 7 is the ordinary polytropic 
index and n is a polytropic index commonly used in as- 
trophysics [58] . The polytropic constant K is given for 
n > d/2 by 



K = 



1 



1 



AS d 2 d / 



,_ 1 r(d/2)r(i - d/2 

r(n+l) 



and for n < — 1 by 
K = - 1 



1 



r( 



(28) 



(29) 



The polytropic constant K plays the role of the temper- 
ature Ti SO in isothermal systems (q = 1, n = 00, 7 = 1) 
and it is sometimes called a "polytropic temperature" . For 
polytropic distributions, the relation between the poly- 
tropic temperature K and the thermodynamical temper- 
ature T = 1/(3 is 



3.4 Other expressions of the distribution function 

We can write the polytropic distribution function (1171) 
in different forms that all have their own interest. This 
will show that different notions of "temperature" exist for 
polytropic distributions. 

(i) Thermodynamical temperature T = 1/(3: the poly- 
tropic distribution ([17)) directly comes from the variational 
principle (fl"6|) . Therefore, (3 = 1/T is the Lagrange mul- 
tiplier associated with the conservation of energy. This is 
the proper thermodynamical temperature to consider, i.e. 
(3 = (dS/ 8E)m is the conjugate of the energy with respect 
to the entropy. Note, however, that T = 1/(3 does not have 
the dimension of an ordinary temperature (squared veloc- 
ity). _ 

(ii) Dimensional temperature 8 = 1/b: we can define a 
quantity that has the dimension of a temperature (squared 
velocity) by setting b = (3/q/j. If we define furthermore 
/» = fi 1 ^ q ~ 1 \ the polytropic distribution ([T7|) can be 
rewritten 



where C n is given for n > d/2 by 



(30) 



/ = /* 



l-b{q-l)t 



(34) 



2(n+ 1)5S=3 
2n - d + 2 



S,2 d ' 



! _ 1 r(d/2)P(l - d/2 



r(n + i) 



2_ Comparing Eq. (J33} with Eqs. (O and ([201) we find that 



e m = l/[6(<7 - 1)] and A = f*/e]l {q 1} for n > d/2 and 

1/(9- 



(31) 



and for n < — 1 by 
2[-(n+l)p 



e = l/[b{\ - q)} and ^ = /*/eo " for n < -1. Sub- 
stituting these expressions in Eqs. ([23)) and ([25)) . we find 
that the relation between the density and the potential 
can be written 



-2n + d- 2 



, d/ 2-i r(d/2)r(-n) 

d r(d/2 - n) 



(32) 



We note that K is a monotonically increasing function of 
T. This remark will have some importance in the follow- 



3.3 Polytropic distributions in physical space 

The density is obtained by integrating Eq. ([17)) on the 
velocity leading to Eqs. j23]) and (T23J). Using Eqs. (f25J) and 
()2l?l) , we find that the density is related to the potential 
<2>(r) by 



1-&(<Z-1)<2> 

where p* is given for n > d/2 by 

S d , (2n-d\ d/2 r{d/2)r{\- d/2 + n) 



(35) 



2 J \ b 



and for n < — 1 by 



_ S± f ( d-2n 
P * ~ ~2 U \~b~ 



d/2 



r{n+l) 



r(d/2)r(-n) 
r{d/2-n) 



(36) 



(37) 



P(r) = 



A - ^*(r) 



(33) 



where A = e m /(K(n+l)) for n > d/2 and A = e /(-if (n+ 
1)) for n < —1. As noted in [37] . the polytropic distribu- 
tion in phase space / = /(e) given by Eq. (|17p has the 
same mathematical form as the polytropic distribution in 
physical space p = p{<P) given by Eq. ([3"3")) with 7 playing 
the role of q and K playing the role of T = 1/(3. In this 
correspondence, 7 is related to q by Eqs. ([2"7)l and lfl"%)) 
leading to 7 = ((d + 2)g — d)/(dq — d + 2) and -fT is related 
to T by Eq. ([30]) . Polytropic distributions (including the 
isothermal one) are apparently the only distributions for 
which /(e) and p(^) have the same mathematical form. 



(iii) Polytropic temperature K: eliminating the poten- 
tial between Eqs. ([TH]) and ([ZD, and between Eqs. (gO)) 
and ([25)) . and using Eqs. (f28| and ([29)). we can express 
the distribution function in terms of the density accord- 
ing to 



1 



p(r) 



l/n 



v 2 /2 



(n + 1)K 



n-d/2 



(38) 



where Z is given for n > d/2 by 



z = Sj2 ^- ^(^w-y+") [g(n + irV39 > 
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and for n < — 1 by 

Z = S d 2 d '^ r{ -^ {d,2 \ -K(n + 1)}^. (40) 
1 {d/2 — n) 

This is the polytropic counterpart of the isothermal dis- 
tribution. The constant K plays a role similar to the tem- 
perature Ti SO in an isothermal distribution. In particular, 
it is uniform in a polytropic distribution as is the temper- 
ature in an isothermal system. This is why K is sometimes 
called a polytropic temperature. 

(iv) Kinetic temperature T k i n (r): the kinetic tempera- 
ture is defined by 



T km (r) = -Av 2 )(r)=p(r)/p(r). 



(41) 



For a polytropic distribution, using the equation of state 
{2Z|), we get 



T kin {r) =Kp(r 



,7-1 



(42) 



For a spatially inhomogeneous polytrope, the kinetic tem- 
perature 7fcj n (r) is position dependent and differs from 
the thermodynamical temperature T = 1//3. The velocity 
of sound c 2 s = p'(p) is given by 



cfr) = K 1P {v) 



7-1 



jT kin (r). 



(43) 



It is also position dependent and differs from the velocity 
dispersion (they differ by a factor 7). Using Eq. (|4"Tj) . the 
distribution function (j3"8")l can be written 



f = B n 



Br, 



p(r) 



[2nT kin (v)}^ 



v 2 /2 



(n+l)T Mn (r) 



-d/2 



(44) 



r(n + 1) 



r(l-d/2 + n)(n + l) d / 2 



(n > d/2) (45) 



r(d/2-n) 

Note that for n > d/2, the maximum velocity can be ex- 
pressed in terms of the kinetic temperature by 



v m (r) = y/2(n+l)T kin (r). 



Using r(z + a)/r{z) ~ z a for z — > +00, we recover the 
isothermal distribution for n — > +00. On the other hand, 
from Eqs. (|4"2"j) and (J33J) , we immediately get T k i n (r) = 
K{\ - (7 - l)<P(r)/K~f) so that 



7- 1 

VT fcm = -i W. 

7 



(48) 



This shows that, for a polytropic distribution, the kinetic 
temperature (velocity dispersion) is a linear function of 



the potentiaQ The coefficient of proportionality is re- 
lated to the polytropic index by (7 — = l/(n + 1) = 
2(q - l)/[(d + 2)q - d]. The relation {3S]) can also be ob- 
tained directly from Eq. (fT^|) [or Eq. ([201) ] noting that 



/ = A(e m - «P) 



i-d/2 



1 



i; 2 /2 



i-d/2 



(49) 



and comparing with Eq. (|44[) . 

(v) Energy excitation temperature T(e): for any dis- 
tribution function of the form / = /(e), one may define 
a local energy dependent excitation temperature by the 
relation 



1 



T(e) 



din/ 
de ■ 



(50) 



For the isothermal distribution, T(e) coincides with the 
usual temperature Ti SO , For the polytropic distribution 
(fP7|) . T(e) = qn/ (3 — (q — l)e. This excitation tempera- 
ture has a constant gradient dT /de = 1 — q related to the 
polytropic index. The other parameter [i is related to the 
value of the energy where the temperature reaches zero. 



3.5 Entropy and free energy in physical space 

Substituting the polytropic distribution function (|38[> in 
the Tsallis entropy (fTTj) . we find after some calculations 
(see Appendix [C]) that 



(fi J 'pdr-f^M^, (51) 



On the other hand, the energy (fUij) can be written 

E=^[pdr+±[ f&dv. (52) 

Therefore, the free energy F = E — TS is given by 

F = i J p<Pdr + n J pdr, (53) 

up to unimportant constant terms. It can also be written 



F =2 



IJ^fr+KJfjfr- pJpV) dr, (54) 



where 



7 _1 _ 7-1 T q _! 

Po " K q-\ h ■ 



(55) 



(47) 

This can be viewed as a free energy F = E—KS associated 



with a Tsallis entropy in physical space S = J (p 1 

pl~ 1 p)dr where 7 plays the role of q and K the role of 
T. Again, it is interesting to note that, for polytropes, the 
free energy in phase space F[f] has a form similar to the 
free energy in position space F[p] [4"7] . 

17 For any distribution of the form / = /(e), we have p = 
p(^) and p = p(^) so that the kinetic temperature (velocity 
dispersion) Tfc in = p/p is a function T^ in — Tfc in (^) of the 
potential. For a polytropic distribution function, this relation 
turns out to be linear. 
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4 Application to the HMF model 

4.1 Generalities 

The HMF model is a system of N particles of unit mass 
m = 1 moving on a circle and interacting via a cosine 
potential. The dynamics of these particles is governed by 
the Hamilton equations |6l8l9j : 



dt 



dH 

dvi ' 



dvj 
dt 



dH 

'Wi' 



H = 



1 



JV 



£- 2 



k 

4tt 



^cos(^-^), 



(56) 



where 6*.; <G [— 7r,7r] and — oo < Vi < oo denote the position 
(angle) and the velocity of particle i and k is the coupling 
constant. The thermodynamic limit corresponds to N — > 
+oo in such a way that the rescaled energy e = 8TrE/kM 2 
remains of order unity (this amounts to taking k ~ 1/N 
and E/N ~ 1). In that limit, the mean field approximation 
becomes exact. The total mass and total energy are given 
by 



M = / pd6, 



E 



f— dOdv 
J 2 



p<Pd6 



W, 



(57) 



(58) 



where f{0,v) is the distribution function, p{9) = J f dv 
the density, E kin the kinetic energy and W the potential 
energy. The potential is related to the density by 



m = - 



k 

2^ 



2tv 



cos(9-9')p(9')d9', 



(59) 



and the mean force acting on a particle located in 9 is 
(F) = —<P'(9). For a cosine interaction, the potential can 
be written 



-B x cos( 



By sin I 



where 



k 

Br ~ 

2tt 



2^ 



B„ 



k 

2^ 



p{9) coa ede, 



p{6) smOdO, 



(60) 



(61) 



(62) 



are the two components of the magnetization vector up to 
a multiplicative factor (note the change of sign with re- 
spect to 0). Substituting Eq. (|6"0"|) in the potential energy 
(j58]) . we find that the energy can be rewritten 



/— dOdv 
J 2 



ttB 2 



(63) 



where B 2 = B 2 + B 2 . It is a nice property of the HMF 
model that the potential energy can be expressed in terms 
of the order parameter B played here by the magneti- 
zation. This is not the case for more complex potentials 
such as the gravitational potential. Note, finally, that the 
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Fig. 1. Density profile for the polytropes n = 1. We have 
represented incomplete (x > x c ) and complete (x < x a ) poly- 
tropes. 



kinetic temperature can be expressed in terms of the pres- 
sure (|2"Tj) so that 



E=l 
2 



27T 



pdO- 



nB 2 



= E k 



kin 



w. 



(64) 



4.2 Complete and incomplete polytropes 



Let us now study polytropic distributions in the context 
of the HMF model. According to Eq. (|33[) . their density 
profile is given by 



p(0) = A 



7-1 



i 

7-1 
+ 



(65) 



where A = A 1 ^ 7 " 1 ). We can assume without loss of gen- 
erality that the distribution is symmetrical with respect 
to the x axis (9 = 0). In that case, the potential can be 
written 

<P = -Bcos9, (66) 

where 



k f 2 * 
B = — p{9) cos 9 d9, 
2tt Jo 



(67) 



is the magnetization (B y = and B = B x ). The density 
profile can be rewritten 



p(9) = A 



where 



7-1 



B 



- KAr -r- m 

For 7 — > 1, we recover the isothermal distribution (20) 
and the relation x = j3B of [5]. For x > 0, the density 
profile is concentrated around 9 = and for x < 0, we get 
a symmetrical density profile concentrated around 9 = ir. 

Like for the gravitational interaction, we can have com- 
plete or incomplete polytropes [JB]. By definition, com- 
plete polytropes have a compact support (the density 
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drops to zero at 9 C < n) while incomplete polytropcs ex- 
tend over the whole domain. Let us first consider the case 
n > 1/2 (i.e. 1 < 7 < 3). The density profile can be 
written 

/ \ 71 

(70) 



p(0)=A 1 + — -cos! 

n + 1 



By symmetry, we can restrict ourselves to the interval < 
9 < 7T and to x > 0. The density profile is monotonically 
decreasing and it has a compact support iff the term in 
parenthesis is negative for 9 — n. This occurs iff 



x > x c = n + 1 



7 



7-1 



(71) 



If x < x c , the density is strictly positive on < 9 < n 
(incomplete polytrope). If x > x c , the density vanishes 
for 9 > 9 C (complete polytrope) such that 



1 



1 



71+ 1 

Therefore, 9 C is given by 



x cos 9 C = 0. 



arccos — 



(73) 



Some typical polytropic profiles are represented in Fig. [TJ 
We now consider the case n < — 1 (i.e. < 7 < 1). The 
density profile can be written 



.4 



P(<>) = — 
and the solution is physical iff 



x x * 



-(n+l) = - 



7-1 



(74) 



(75) 



In that case, the polytropes are always incomplete. 



4.3 The magnetization and the polytropic temperature 

For the density profile the mass and the magnetiza- 
tion are given by 

M = A [ (l + ^— UcosflV 1 d6, (76) 



kA r ( -7-1 \~ 

B= — / 1+- xcos9) cos9d9. (77) 

2 Wo \ 7 J + 

It is useful to introduce the 7-deformed Bessel functions 

1 f 27T ( 7 - 1 \^ 
I-y mix) = — / 1 H cccos^ cos(m9) d9 . 

^ Jo \ 1 J + 

(78) 

For 7 — > 1, we recover the Bessel functions I m {x). For 
l<7<3(n>l/2) and x < x c , we have 

1 r ( 7- 1 \ T ^ T 

I-y m(x) = — / f H x cos 9 cos(m9) d9, 

2tt J-k V 7 / 

(79) 



-o 




Fig. 2. Magnetization as a function of the polytropic temper- 
ature for n > 1/2. 



(72) and for x > x c , we have 



1 f ( 7-1 \~ 

I-ym(x) = — / 1H x cos 9 J cos(m9)d9. 

2tt 7_e c \ 7 / 

(80) 

We also have I 1>m {— x) = I Jtm (x) if m is even and 
Xy )TO (— x) = —I-y, m (x) if m is odd. In terms of these in- 
tegrals, we can rewrite the relations ([75|) and ([77)1 in the 
form 

A=-^-, (81) 



and 



27r/ 7i0 (a;) 
2ttB _ / 7 ,i(x) 



(82) 



kM I-y,o(x) 

Combining Eqs. ((5T|) and ([52^1 and introducing the 

normalized polytropic temperature [9]: 



V = (2tt) 



471"^ : 



we find that 



2 7 7fl (x) 



(83) 



(84) 



Equations (|82|) and (|84|) determine the magnetization i? 
as a function of the polytropic temperature K in a para- 
metric form (with parameter x). For 7 — > 1, we recover 
the self-consistency relation (26) of [5]. Some curves are 
plotted in Fig. [5] for n > 1/2. The case of negative indices 
n < — 1 is very similar to the isothermal case [H] and will 
not be illustrated in detail. For x — s- x*, b — > 1, r/ — > +00 
and e — > —2. 



4.4 The energy 

4.4.1 The homogeneous phase 

In the homogeneous phase, the density is uniform with 
value p = M/ (2tt) and the magnetization vanishes (B = 
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Fig. 3. Magnetization as a function of the energy for n > 1/2. 



0). Therefore, the energy reduces to the kinetic energy 

1 f 2n 1 f 27T f M\~< 

E = -J o P d6=- 2 Kj o P^0 = K^-) . (85) 

Introducing the dimcnsionlcss energy [5]: 

8ttE 



kM 2 ' 



(86) 



and recalling the expression (j8"3")l of the dimensionless 
polytropic temperature, we find that the caloric curvcP^I 
for the homogeneous distribution is simply [9]: 



V = 



(87) 



It has been shown in 9 48] that the homogeneous phase 
is stable iff 



V < r lc = l, e > e c = 



4.4.2 The inhomogeneous phase 



1 



To compute the kinetic energy in the inhomogeneous 
phase, without introducing new integrals, we shall use a 
trick inspired from the astrophysical problem [46j . The 
condition of equilibrium can be written as a condition of 
hydrostatic balance (see Appendix [5} [3] : 



dp _ d<£> 
d0 = ~ P ~dl' 



(89) 



For a polytropic equation of state, we have 

p = Kpf. (90) 
Therefore, the foregoing equation can be integrated into 

<P = — ^-p^- 1 + C. (91) 
7-1 



18 We shall explain in Sec. [5] why this is the relevant caloric 
curve to consider. 



We must be careful that this relation is valid only in the 
region where p > 0. We must therefore consider two cases. 
For incomplete polytropes, using <P(ir) = B, we get 



-(py-i-p^y-^ + B. 



(92) 



On the other hand, for complete polytropes, using p(0 c ) 
and $(6 C ) = -Bcos0 c , we get 



<P = — I —p~<- 1 -Bcos9 c 
7- 1 



(93) 



These two cases can be combined in a single formula 



<2> = — 



7-1 



{p 1 -' -pin) 1 -') -Boost 



(94) 



For incomplete polytropes 8 C = n, and for complete poly- 
tropes p(ir) = 0. Substituting Eq. (|M|) in the potential 
energy ([55)1 and using Eq. (f5D|) . we obtain 



W = - 



1 7 



2 7 -l 



P d0+ l-^Mpiny- 1 



27-I 
1 



MB cos 9 r 



(95) 



In the first integral, we recognize the kinetic energy (|6~H) . 
Therefore, we obtain 



ki7i 



7-1 



w 



1 



KMp{ir) 



7-1 



1 7 - 1 



MB cos 6 r 



(96) 

Using the expression (jMJ) of the potential energy, we get 

'V — 1 7rR 2 1 'Y — 1 1 

E km = -i — - MB cos0 c + -KMp{-Ky~\ 

7 2 7 2 

(97) 



(88) Finally, the total energy E = E kin + W is 

£ = - l-^l M Bcose c + i|fMp(7r)' ir - 1 . (98) 

7/c 27 2 

Using Eqs. (|82p and (|6"B"|) . the normalized energy (|5rJ|) can 
be written 



e = --b 2 - 2^ — -bcos8 c + 2- f 1 - - — -ac 

7 7 x \ 7 



(99) 



Using Eq. (1821) . it can also be written explicitly in terms 
of x as 

2 1^ 1 {x) 2 j-l /^) 

— Z ; — 7 COS U r 



1 1~f,oi x ) 2 7 Aso(z) 

, 2/ 7 ,i(ar) 



a; 7 7 ,o(x) 



' 1 - 



7 



(100) 



Equations ([52^1 and poop determine the magnetization 6 
as a function of the energy e in a parametric form. Some 
curves are plotted in Fig.[3]for n > 1/2. On the other hand, 
the caloric curve 77(e) is determined by Eqs. (|84p and (|100j) . 
Some curves are plotted in Fig. [4] for n > 1/2 and in Fig. 
[5]for n < —1. They complete Fig. 12 of [3] by adding the 
inhomogeneous branches. They will be analyzed in Sec. [5] 
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Fig. 4. Caloric curve for polytropes with positive indices n > Fig. 6. Inverse averaged kinetic temperature as a function of 
1/2. the energy for polytropes with positive indices n > 1/2. 
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r|=l/e 
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and use Eqs. and (1521 . wc obtain 
w = -2b 2 . 



(106) 



Finally, the total energy E = Eki n + W can be written in 
dimensionlcss form 



e = (-kin + w = eicin - 2b 2 . 
From Eqs. j99| and (|T07| . we obtain 



(107) 



l^±b 2 2 1 —! 



bcost 



Fig. 5. Caloric curve for polytropes with negative indices n < 
-1. 



4.5 The averaged kinetic temperature 

The local kinetic temperature is defined by 

/ fv 2 dv 



7 7 

+2 b - (l - l=± x 
x V 7 



(108) 



Using Eq. (j52"j) , this can also be written explicitly in terms 
of x as 

e km =0 = 2^^S - 2 7^1^iM cosr 



T km (0) = (v 2 )(6) = 
The averaged kinetic temperature is therefore 
(Tun) = ± P (0)(v 2 )(8)d6= 2 



AI 



AI 



(101) 



(102) 



7 I-y,o( x ) 2 7 I-yA x ) 
+ 2U)f, 7-1 
x I 7j o(a;) 



7 ; + 



(109) 



Defining the normalized kinetic energy and the normalized 
averaged kinetic temperature by 



8-KEki 







47r(T, 



A- 



fcM 2 ' fcM 
we obtain 

= Efein- 

If we define the normalized potential energy by 

8nW 

w 



kM 2 ' 



(103) 



(104) 



(105) 



Equations (|109p and (|100p determine the physical caloric 
curve efci n (e) or 0(e) in a parametric form. Some curves 
are represented in Fig. |5]for n > 1/2. In the homogeneous 
phase, we have 

e = e fem = 0. (110) 
These physical caloric curves will be analyzed in Sec. El 



4.6 The free energy 

According to Eq. ([55|. the free energy of the polytropes 
can be written in terms of the density as 



F = - [ p<Pdd + -^— 



P 1 dQ, (111) 
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up to an unimportant constant. In the second term, we 
recognize the kinetic energy (|64|) while the first term is 
just the potential energy, so that 



F = W 



-Eki 



7-1 

Introducing the normalized free energy 

8ttF 



we obtain 



/ 



f = w + 



kM 2 ' 



7-1 



(112) 



(113) 



(114) 



Substituting Eqs. ([TPS)) and ([TU5| in Eq. ([TUij) . we find 
that 



/ 



2(2- 7) b 2 



-b cos 9 C 



' ' \ - ————x 



7 



j — 1 x 



(115) 



Using Eq. 

x as 



this can be explicitly written in terms of 



2(2-7) J 7;1 (x) 2 4 J 7 ,i(a;) c| 
7 I 7 ,o(x) 2 7 Iy,o{x) 
1 4 J 7 ,i(i) ( x _l_J_ x 



7 - 1 x I lfi {x) 



7 



(116) 



Equations (|116[) and determine the free energy as a 
function of the polytropic temperature f(rj) in a paramet- 
ric form. Some curves are represented in Figs. [7] and [5] for 
n > 1/2. They will be analyzed in Sec.[5l In the homoge- 
neous phase, Eq. (|112p reduces to F = 2E/ (7 — 1). Using 
Eq. dHZJ) , wc get 



/ = 



2 1 



7- in 



(117) 



0.65 - 




0.55 - 



Fig. 7. Free energy as a function of the inverse polytropic 
temperature for n = 0.7. It exhibits a first order canonical 
phase transition marked by the discontinuity of df / dr\ at r\ — 
Vt- 



<*-2- 




Fig. 8. Free energy as a function of the inverse polytropic tem- 
perature for n — 2. It exhibits a second order phase transition 
marked by the discontinuity of d 2 f /drj 2 at 77 = rj c . 



Defining the normalized entropy by 



4.7 The entropy 

According to Eq. (f51"j) , the entropy of the polytropes can 
be written in terms of the density as 

S=-(n-^pJpd6, (118) 

up to an unimportant constant. In the second term, we 
recognize the kinetic energy so that 

S = -2 - pE kin . (119) 

On the other hand, the inverse thermodynamical temper- 
ature (3 = 1 jT is related to the polytropic temperature K 
by Eq. (|30p . Therefore, the entropy can be rewritten 

s ~H)^- (120) 



A/— 

we obtain 

s = - (n - ^ e km T]^. (122) 

Using Eqs. ([51)1 and (|108p . this gives s as a function of 
x. Finally, Eqs. (|122[) and (|100[) determine the entropy as 
a function of the energy s(e) in a parametric form. Some 
curves arc represented in Figs. |pj and [TUJ for n > 1/2. They 
will be analyzed in Sec. [SJ For the homogeneous phase, 
using e = efc !rl and r\ — 1/e, we get 




(123) 
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Fig. 9. Entropy as a function of energy for n = 0.53. It exhibits 
a microcanonical first order phase transition marked by the 
discontinuity of ds/de at e = et. To avoid the occurrence of 
large numbers when n ~ 0.5, we have plotted s* = — c fein 
and s » omo = — 1/e as a function of e. 
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Fig. 10. Entropy as a function of energy for n — 0.7. It ex- 
hibits a microcanonical second order phase transition marked 
by the discontinuity of d 2 s/de 2 at e = e c . There is also a small 
convex dip (hardly visible) evidencing a region of negative spe- 
cific heats and a canonical first order phase transition. 



5 The thermodynamical caloric curve 
5.1 General comments 

Basically, we want to study the optimization problems Q 
and §5§ where S is the Tsallis entropy (fTTj) . The maxi- 
mization problem ((4]) corresponds to a condition of micro- 
canonical stability and the minimization problem (JS|) to a 
condition of canonical stabilit\P^I. It is well-known that the 
two optimization problems ((4]) and ([5]) have the same crit- 
ical points canceling the first order variations of the ther- 
modynamical potential. The critical points of the maxi- 
mization problem (j4]) are given by Eq. (|16p where (3 and 



19 Recall that these optimization problems can be given a dy- 
namical interpretation (see Sec. 12. 2\ . For convenience, we use 
here a thermodynamical language based on the thermodynam- 
ical analogy described in Sec. 12.21 We should add the prefix 
"pseudo" (i.e. pseudo entropy, pseudo temperature...) but for 
simplicity, we avoid it. 



a are Lagrange multipliers associated with the constraints 
E and M. Here, /3 = 1/T represents the inverse thermo- 
dynamical temperature. It satisfies the thermodynamical 
relation (3 = (dS/ 8E)m ■ For the Tsallis entropy, the criti- 
cal points in the microcanonical ensemble are given by Eq. 
(fTT)) . The critical points of the minimization problem (J5J) 
are given by 5F + aTSM = where a is a Lagrange mul- 
tiplier associated with the constraint M. Obviously, these 
first order variations are equivalent to Eq. ([TB"]) so that the 
critical points in the canonical ensemble are again given 
by Eq. (|T7|) . Although the critical points are the same, 
their stability can differ in the microcanonical and canon- 
ical ensembles. From general theorems [50], we know that 
a critical point is stable in the microcanonical ensemble 
(maximum of S at fixed E and M) if it is stable in the 
canonical ensemble (minimum of F at fixed M) but the 
reciprocal may be wrong. In that case, we have ensembles 
inequivalence. Related to this property, we can show that 
the thermodynamical specific heat C = dE/dT is always 
positive in the canonical ensemble while is can be posi- 
tive or negative in the microcanonical ensemble. These are 
general results valid for all entropies of the form ©, in- 
cluding the Tsallis entropy , provided that we use the 
proper thermodynamical variable E and (3 = (dS/dE) m 
that are canonically conjugate with respect to the entropy 
S. Therefore, the proper thermodynamical caloric curve 
associated with the optimization problems (j4|) and ([5]) is 
the one giving the inverse temperature (3 = 1 /T (Lagrange 
multiplier) as a function of the energy E. It is from this 
curve that we can study ensembles inequivalence [5UJ. It 
is also from this curve that we can apply the Poincare 
theory on the linear series of equilibria [59160] . Now, for 
polytropic distributions, we have seen in Eq. ([50]) that the 
polytropic temperature K is always a monotonic func- 
tion of the thermodynamical temperature T. Therefore, 
for convenience, we can cquivalently consider the caloric 
curve giving the polytropic temperature K as a function 
of the energy E (the preceding properties are unchanged) . 
This is the choice that has been made in the gravitational 
problem [46139] and that we shall make here. Therefore, 
we shall call K(E) the thermodynamical caloric curve. 
In dimensionless variables, it corresponds to 77(e). Several 
caloric curves 77(e) have been plotted in Figs. [4] and [5] for 
different indices n. In order to understand their behav- 
ior let us first consider the asymptotic limits x — > and 



5.2 The limit x 

For x — > 0, the deformed Bessel functions (]78]) can be 
approximated by 



I-t,o( x ) = 1 + 



2-7 2 
4 7 2 X 



(124) 



(125) 
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Then, we obtain the expansions 



5.3 The limit x 



-co 



x ( 7 -2)(2 7 -l) 3 
2 7 + 16 7 3 



2-7 2 
77 = 7 H — x 



1 2 - -35" 



£=- + (27 -57-2)-^ + 

7 87 J 

1 o X 2 



Efcin = @ = 



7 



(2 7 2 - 7 - 2) 



8 7 3 



(126) 
(127) 
(128) 
(129) 



From these expressions, we can draw the following conclu- 
sions. Let us introduce the critical indices 



1 



2.8507811. 



0.54031242. 



(130) 



(131) 



Let us first consider the case n > 1/2 so that 1 < 7 < 3. 
The quantity 27 2 — 57 — 2 is positive for 7, < 7 < 3, 
i.e. 1/2 < n < n*, showing that the energy increases as 
x increases. The quantity 27 2 — 57 — 2 is negative for 
1 < 7 < 7* j i-e. n > n* . In that case, the energy decreases 
as x increases. On the other hand, for 2 < 7 < 3, i.e. 
1/2 < n < 1, the inverse temperature 77 decreases as x 
increases. By contrast, fori < 7 < 2, i.e. n > 1, the inverse 
temperature rj increases as x increases. This explains the 
behavior of the caloric curve 77(e) close to the bifurcation 
point. Close to that point, we have 



V - 7]c 



(2-7)7 2 
2 7 2 - 57 - 2 



The specific heat is 



dE 2 dE 



(132) 



(133) 



Since T is a monotonically increasing function of K , we 
can equivalently consider the quantity 



c K = -n 



de 
drj 



At the bifurcation point, we have 



C K = 



2 7 2 - 57 - 

2^7 



(134) 



(135) 



The specific heat close to the bifurcation point is positive 
for 1/2 < n < n*, negative for n* < n < 1 and positive 
again for n > 1. It vanishes for n = n* and it is infinite 
for n=l. Let us now consider the case n < — 1 so that 
< 7 < 1. In that case, 77 increases as x increases. On the 
other hand, the quantity 27 2 — 57 — 2 is always negative 
so that e decreases as x increases. Therefore, the specific 
heat is always positive. 



Let us first consider n > 1/2. For x — > +00, we find that 
9 C — > ^. The deformed Bcssel functions can be approxi- 
mated by 



/ 7,o( X ) 



1 



1 



2tt (n + 1)" 
1 1 



27r(n + l)« v r(2±™)' 
The normalization constant (I5H behaves like 



(136) 
(137) 



1 1 r(^) 

A ~ M (n + 1) " — ^= 2 



and the density profile 
p{6) = 



(138) 

tends to the limit distribution 
(139) 



Finally, the thcrmodynamical parameters tend to the limit 
values: 



1 r 



. 2 > 



2 1 ' (v^) 1 /" r(2±H)^ 



e ->■ £00 = -- 



r 



272 



r 



r( 2± 2; 



Chi 



e 



(140) 
(141) 

(142) 
(143) 



This implies in particular that the polytropes n > 1/2 
only exist for e > and 77 < 7^00. Beyond these points, 
there is no polytropic distribution anymore. In that case, 
the system is likely to converge towards the Lynden-Bell 
distribution (see Sec. 17. 2|) . 



5.4 Phase transitions 

We are now in a position to describe in detail the phase 
transitions that occur in the thermodynamical caloric 
curve K(E) or, in dimensionless form 77(e), plotted in Figs. 
S] and [SJ Note that this curve represents the full series of 
equilibria containing all critical points of entropy or free 
energy. We now have to determine which solutions corre- 
spond to fully stable (S) solutions (global maxima of S 
in MCE or global minima of F in CE), metastable (M) 
solutions (local maxima of S in MCE or local minima of 
F in CE) or unstable (U) solutions (saddle points of S or 
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Fig. 11. For n > uctp = 1 (specifically n = 2), the 
caloric curve exhibits second order canonical and microcanon- 
ical phase transitions. The ensembles are equivalent. 



F). To make this selection, we can compare the equilib- 
rium entropy or free energy of the solutions in competition 
(see Sees. 14.61 and !4.7p or use the Poincare theory of linear 
series of equilibria |59l60j . 

Let us first consider the canonical ensemble (CE) 
where the control parameter is the inverse polytropic tem- 
perature 77. For n > tictp = 1, the caloric curve e{rj) ex- 
hibits a second order phase transition (see Figs. 21 [5] and 
Hip marked by the discontinuity of de/drj at 77 = r\ c . The 
homogeneous solution is fully stable for rj < rj c and it be- 
comes unstable for 77 > r\ c . At that point it is replaced by 
the inhomogeneous branch that is fully stable. The series 
of equilibria presents one solution for 77 < r\ c (the solution 
6 = that is fully stable) and three solutions for 77 > ?7 C 
(two solutions ±6 7^ that are fully stable and one solu- 
tion 6 — that is unstable); see Fig. [2] For n = uctp = lj 
there is a degeneracj0 at 77 = ?; c = 2 since there exists an 
infinity of solutions — 1/2 < 6 < 1/2 with the same free 
energy (see Figs. [21131 and 112)1 ■ For 1/2 < n < ii C tp = 1, 
the caloric curve e(rj) exhibits a first order phase transi- 
tion between the homogeneous and inhomogeneous phases 
(see Figs. 0] [7] and [T5|) marked by the discontinuity of 
e at 77 = rj t and the existence of metastable branches. 
The series of equilibria presents one solution for 77 < 77* 
(the solution b = that is fully stable), five solutions for 
77* < ?7 < ?7 C (two solutions ±6 7^ that are fully stable or 
metastable, two solutions ±6' 7^ that are unstable and 
one solution 6 = that is fully stable or metastable) and 
three solutions for r/ > rj c (two solutions ±6 7^ that are 
fully stable and one solution 6 = that is unstable); see 
Fig. [21 Therefore, tictp — 1 and rj c = 2 (corresponding 
to e c = 1/2) is a canonical tricritical point separating first 
and second order phase transitions. The canonical phase 
diagram in the (n, 77) plane is represented in Fig. 1141 

Let us now consider the microcanonical ensemble 
(MCE) where the control parameter is the energy e. For 
n > umcp — 0.68, the caloric curve 77(e) exhibits a sec- 
ond order phase transition marked by the discontinuity of 



This is a bit similar to critical polytropes n c = 3 in astro- 
physics [61] . 




Fig. 12. For n = uctp = 1, there is a degeneracy in the canon- 
ical ensemble for 77 = r\ c = 2. In the microcanonical ensemble, 
there is a region of infinite specific heats. 



2.6 


_ 1 1 

n=0.7 




\ 




2.5 






U \ 




2.4 


^\ S 

\. E i 




u , 


\ \ 

\ M - 
E 2 \ 




\ ^vT 








2.3 


M 
1 




' u 


s y 



\ I I I I I I 

0.25 0.3 0.35 0.4 

e 



Fig. 13. Caloric curve in the canonical ensemble for 0.5 < n < 
tictp = 1 (specifically n = 0.7). It exhibits a canonical first or- 
der phase transition marked by the discontinuity of e at 77 = 774. 
The temperature of transition can be obtained by comparing 
the free energy of the homogeneous and inhomogeneous phases 
(see Fig. [7|. The Maxwell construction does not directly apply 
in the present situation since we have used 1/K instead of j3 in 
the caloric curve. The stability of the solutions (denoted S, M, 
and U) can be settled by using the Poincare theorem [59 60 . 
The temperatures 77* and r\ c are spinodal points marking the 
end of the metastable branches. 



dri/de at e = e c (see Figs. H MUD 03 and E3). The homo- 
geneous solution is fully stable for e > e c and it becomes 
unstable for e < e c . At that point it is replaced by the inho- 
mogeneous branch that is fully stable. For n > tictp = L 
the specific heat is positive and the ensembles are equiv- 
alent (see Fig. QT]). For n^icp < n < uctp = 1, there 
is a region of negative specific heat in the microcanoni- 
cal ensemble between e' and e c (see Fig. IT5j) that is re- 
placed by a phase transition in the canonical ensemble 
(see Fig. [T2]). The series of equilibria presents one solu- 
tion for e > e c (the solution 6 = that is fully stable) 
and three solutions for e < e c (two solutions ±6 7^ that 
are fully stable and one solution 6 = that is unstable); 
see Fig. [3] The corresponding microcanonical phase di- 
agram in the [n, e) plane is represented in Fig. 1161 For 
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0.5 0.6 0.7 0.8 0.9 1 1.1 12 1.3 1.4 1.5 
n 

Fig. 14. Phase diagram in the canonical ensemble. There is 
a canonical tricritical point located at ncrp — 1 and r\ — 2 
separating first and second order phase transitions. 



timtp — 0.563 < n < tlmcp — 0.68, there is a very in- 
teresting situational in which the caloric curve exhibits a 
second order phase transition at e c between the homoge- 
neous phase and the inhomogeneous phase (as before) and 
a first order phase transition at et between two inhomoge- 
neous phases (see Fig. [17]). The first order phase transition 
is marked by the discontinuity of r\ at e = t t and the ex- 
istence of metastable branches. For n = 0.6, the series of 
equilibria presents one solution for e > e c (the solution 
6 = that is fully stable), three solutions for < e < e c 
(two solutions ±6 ^ that are fully stable and one solu- 
tion 6 = that is unstable), seven solutions for < e < e* 
(four solutions ±6 ^ and ±6" ^ that are fully stable or 
metastable, two solutions ±6' ^ that are unstable and 
one solution 6 = that is unstable) and three solutions for 
e < e'„ (two solutions ±6^0 that are fully stable and one 
solution 6 = that is unstable). Therefore, timcp — 0.68 
and e ~ 0.3492 (corresponding to ?y = 2.336) is a micro- 
canonical critical point marking the appearance of the first 
order phase transition. At n = umtp — 0.563, the ener- 
gies of the first and second order phase transitions coincide 
(e t = e c ) and, for 1/2 < n < umtp — 0.563, there is only a 
first order phase transition at e = et between homogeneous 
and inhomogeneous states. Therefore, umtp — 0.563 and 
e ~ 0.3603 (corresponding to r\ = 2.776) is a microcanoni- 
cal tricritical point separating first and second order phase 
transitions. For rt* ~ 0.54 < n < umtp — 0.563, there re- 
mains a sort of second order phase transition at e = e c for 
the metastable states (see Fig. [T5|l . Between e' + and e c , the 
specific heat is negative. For n = n* ~ 0.54, the energies e'„ 
and e c coincide and for 1/2 < n < n„, the "metastable" 
second order phase transition disappears. In that case, 
there is only a first order phase transition (see Figs. [H] and 
I19|) . The scries of equilibria presents one solution for e > e* 
(the solution 6 = that is fully stable), five solutions for 
e c < e < (two solutions ±6^0 that are fully stable 
or metastable, two solutions ±6' ^ that are unstable 
and one solution 6 = that is fully stable or metastable) 
and three solutions for e < e c (two solutions ±6 ^ that 



21 To our knowledge, this is the first example of that kind. 
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Fig. 15. Caloric curve in the microcanonical ensemble for 
riMCP < n < ncrp ~ 1 (specifically n = 0.7). The inho- 
mogeneous states are fully stable. There is a region of negative 
specific heats between e et e c . In the canonical ensemble, this 
region is replaced by a phase transition. 
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Fig. 16. Phase diagram in the microcanonical ensemble. In 
this range of indices, there is only a second order phase tran- 
sition at energy e c . The states between e et e c have negative 
specific heats. They are inaccessible in the canonical ensemble 
(saddle points of free energy). If we consider only fully stable 
states, the ensembles are inequivalent between ei and £2 (see 
Fig. I13|l . If we take into account canonical metastable states, 
the ensembles are inequivalent between e and e c . 



are fully stable and one solution 6 = that is unstable); 
see Fig. [21 The microcanonical phase diagram in the (n, e) 
plane, summarizing all these results, is represented in Fig. 

M 

In conclusion, for 1/2 < n < tl^itp, we have canon- 
ical and microcanonical first order phase transitions, for 
umtp < n < tlmcp , we have canonical first order tran- 
sitions, microcanonical first order transitions and micro- 
canonical second order transitions, for hmcp < n < 
ncrp = 1, we have canonical first order phase transitions 
and microcanonical second order phase transitions and for 
n > ncrp = 1, we have canonical and microcanonical sec- 
ond order phase transitions. The ensembles are equivalent 
for n > ncrp = 1 and inequivalent for n < ncrp — 1 in 
some range of energies. The microcanonical and canoni- 
cal tricritical points do not coincide. Similar observations 
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Fig. 17. Caloric curve in the microcanonical ensemble for 
timtp < n < riMCP (specifically n = 0.6). Interestingly, there 
exists microcanonical second order and first order phase tran- 
sitions on the same curve. The second order phase transition is 
marked by the discontinuity of dr/fde at e = e c . The first order 
phase transition marked by the discontinuity of n at e = et. 
The stability of the solutions (denoted S, M, and U) can be 
settled by using the Poincare theorem |59I60| . The energies e'„ 
and e» are spinodal points marking the end of the metastable 
branches. 
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Fig. 18. Caloric curve in the microcanonical ensemble for n» < 
n < timtp (specifically n = 0.55). If we only consider fully 
stable states, there is only a first order phase transition at 
e = et. If we consider the metastable states, there is a second 
order phase transition at e — e c . The spinodal points marking 
the end of the metastable branches are e'„ and e* . 
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Fig. 19. Caloric curve in the microcanonical ensemble for 
1/2 < n < n» (specifically n = 0.53). There is only a first 
order phase transition at e = et- The spinodal points mark- 
ing the end of the metastable branches are e c and e,. Note 
that there exists a small region of negative specific heats be- 
tween the turning points of energy and temperature that is 
metastable. 
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Fig. 20. Phase diagram in the microcanonical ensemble. For 
n > riMCP {timcp — 0.68, emcp = 0.3492), there is only 
a second order phase transition. For timtp < n < timcp 
(timtp — 0.563, smtp = 0.36032) there is a second or- 
der phase transition and a first order phase transition. For 
1/2 < n < riMTP, there is only a first order phase transition. 
For n* < n < timtp (n» — 0.54, e = 0.3508), the metastable 
branch presents a second order phase transition. 



6 The physical caloric curve 



have been made for other models [62-65]. For n < — 1, 
we have canonical and microcanonical second order phase 
transitions and the ensembles are equivalent. 

Using the thermodynamical analogy of Sec. 12.21 we 
conclude that the polytropes denoted (S) or (M) that are 
maxima of pseudo entropy S at fixed mass and energy 
are nonlinearly dynamically stable according to the "mi- 
crocanonical" criterion We cannot conclude that the 
polytropes denoted (U) are dynamically unstable since (HJ 
provides just a sufficient condition of nonlinear dynamical 
stability. 



We have explained previously that the proper thermody- 
namical caloric curve to consider when we study phase 
transitions and ensembles inequivalence is the curve T(E) 
where T = is the thermodynamical temperature with 
(3 = {dS/dE)t l [. However, in practice, the caloric curve 
that is experimentally or numerically accessible is the 
curve (Tkin){E) that gives the averaged kinetic temper- 
ature as a function of the energy E. This is the one, for 
example, that is represented in Fig. 4 of Antoni & Ruffo 
[5]. It is therefore useful to study this curve specifically. 
We will call it the physical caloric curve. In dimensionless 
form, it is given by 0(e). Several caloric curves l/(9(e) 
have been plotted in Fig. [S] for different indices n > 1/2. 
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In order to understand their behavior, let us consider the 
asymptotic limit x — > 0. 

From expressions (|128[) and (|129|) . we can draw the 
following conclusions. Let us introduce the critical indices 



7o 



n 



l + s/rf 



17-3 



1.2807764. 



3.5615528. 



(144) 



(145) 



Let us first consider the case n > 1/2 so that 1 < 7 < 3. 
The quantity 2~f 2 — 57 — 2 is positive for 7* < 7 < 3, 
i.e. 1/2 < n < n„, showing that the energy increases as 
x increases. The quantity 27 2 — 57 — 2 is negative for 
1 < 7 < 7„, i.e. n > n», showing that the energy decreases 
as x increases. On the other hand, the quantity 27 2 — 7 — 2 
is positive for 70 < 7 < 3, i.e. 1/2 < n < rig, showing that 
the inverse kinetic temperature decreases as x increases. 
The quantity 27 2 — 7— 2 is negative for 1 < 7 < 70, i.e. n > 
no, showing that the inverse kinetic temperature increases 
as x increases. This explains the behavior of the physical 
caloric curve l/0(e) close to the bifurcation point. Close 
to that point, we have 



o-e c 



2 7 2 - 7 - 2 



2 7 2 - 57 - 
The physical specific heat is 

dE 



kit 1 



M de 

Tdd' 



d(Tkin) 

At the bifurcation point, we have 

M 2 7 2 - 57 - 2 



2 2 7 2 - 7 - 2 



(146) 



(147) 



(148) 



The physical specific heat close to the bifurcation point is 
positive for 1/2 < n < n„, negative for n* < n < no and 
positive again for n > tiq. It vanishes for n = n* and is 
infinite for n = uq. Let us now consider the case n < — 1 
so that < 7 < 1. In that case, the term 27 2 — 7 — 2 is 
always negative so that 1/0 increases as x increases. On 
the other hand, the term 27 2 — 57 — 2 is always negative 
so that e decreases as x increases. Therefore, the physical 
specific heat is always positive. 

Remark: We must be careful that we cannot deduce 
any stability result from this curve. For example, the 
Poincare theorem [59160] and the general results on en- 
sembles inequivalence [50] apply to the thermodynami- 
cal caloric curve T(E), not to the physical caloric curve 
(Tki n )(E). Furthermore, the notion of canonical ensemble 
is only defined in terms of the variables {E,T) via the 
Legendre transform F = E — TS, not in terms of the vari- 
ables (E, (Tkin))- Note, in particular, that for 1 < n < n , 
the solutions close to the bifurcation point have negative 
physical specific heat Chin = dE / d(Tki n ) although they 
are stable both in canonical and microcanonical ensem- 
bles (they have positive specific heat C = dE/ dT Pi 



22 A similar result has been observed recently for the Lynden- 
Bell distribution (in preparation). 



7 The critical index n c = 1 
7.1 Analytical expressions 

We have seen that the index n c = 1 (i.e. j c = 2, q c = 3) 
is particular because it corresponds to a canonical tricriti- 
cal point. Furthermore, we will see that for this particular 
index, the algebra greatly simplifies and analytical results 
can be obtained. This is because the relationship (|6"5"]) be- 
tween the density and the potential is lincail 2 ^!. For n c = 1 , 
the density profile is 



p(9) = A ( 1 + - cos I 



(149) 



For x < x c = 2, the deformed Bessel functions take the 
simple form 



h,o( x ) = !> hAx) = -■ 



(150) 



Then, we get 



M T 

A=—, v = Vc = 2, b=-, (151) 

Z7T 4 



1 X „ 1 X 2 

e= , 9 = - H . 

2 16 2 16 



(152) 



We note in particular that the polytropic temperature is 
constant ?y = n c = 2 (i.e. K c = fc/4 in dimensional form) so 
that there exists an infinity of solutions (the range of mag- 
netizations — 1/2 < b < 1/2) with the same temperature. 
However, they have different energies (see Fig. 2]) ranging 
from 1/4 to 1/2. We will show analytically in Sec. [5] that 
these solutions are stable in the microcanonical ensemble 
while they are degenerate in the canonical ensemble. In the 
range 1/4 < e < 1/2, the thcrmodynamical caloric curve 
7] = t] c has an infinite specific heat C = dE / dT = 00. The 
corresponding physical caloric curve is given by 



0= 1 -e. 



It has a constant negative specific heat 



Chi 



M de 
Yd0 



(153) 



(154) 



23 Note that a similar linear relationship occurs in 2D turbu- 
lence between the vorticity and the stream function for min- 
imum enstrophy states. Now, the enstrophy r-i — j ui 2 dr, 
which is a quadratic functional, can be interpreted as a par- 
ticular Tsallis functional with q — 2 (in our notations). This 
precisely corresponds to the case 7 = 2 considered here (in 2D 
turbulence, we directly work in physical space so that q plays 
the role of 7) . We may note that an approximate linear uj — %p 
relationship has been observed in a plasma experiment [28 29 
and interpreted as a case of incomplete violent relaxation [30] . 
It is interesting to note that the same index 7 = 2 arises in 
the HMF model (see Sec. I7.2[l . Note, however, that in many 
other situations of 2D turbulence, the oj — ip relationship is not 
linear, so that this relation is not universal and may lead to 
physical inconsistencies [30] . 
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This is an example where the thermodynamical caloric 
curve and the physical caloric curve give very different 
results. 

For x > x c = 2, we have cos(9 c = — 2/x and the de- 
formed Bcsscl functions take the form 



h,o( x ) 



1 

2^ 



\/ x 2 — 4 + 2 arccos ( — 



1 

2^ 

Then, wc obtain 



A = 



M 



s/a 



b = 



2 arccos I 



x I 2 
z — 4 + — arccos 



V = 



2/ 2 ,i (x)' 



h,o{ x ) ' 



hAx) 2 | 2 I 2A {x) 
h : o(x) 2 xl 2 fi(x)' 



(155) 
(156) 

(157) 
(158) 
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Fig. 21. Kinetic temperature T as a function of the energy 
U for n = 1. Our theory of polytropes seems to explain the 
anomalies displayed on the numerical caloric curve of Antoni 
& Ruffo [8] close to the critical energy, unlike the Lynden-Bell 
approach which leads to isothermal distributions. 



n hAx) 2 2 J 2 , 1 (aQ 

= 7 , , 2 + ~ T 1 v (159) 

For x — > +00, we have the equivalents l2,o{x) ~ an d 
-^2,1 (x) ~ |. Therefore, the normalization constant be- 
haves like A ~ — and the density tends to the limit dis- 
tribution 

0(0) = yCOB0, (160) 

for 9 < 7r/2. On the other hand, the thermodynamical 
parameters tend to 

•?->4, e ->-^, 0^. (161) 

7.2 Application to the case Mq = 1 

We shall now apply our theory of polytropes to interpret 
the physical caloric curve shown in Fig. 4 of Antoni & 
Ruffo [5] . It is obtained from an initial waterbag distribu- 
tion with magnetization Mo = 1 and it has been confirmed 
by other groups [13114) . Up to date, there is no explana- 
tion of the "anomalies" that take place close to the critical 
energy U c = 3/4. These anomalies are: (i) the bifurca- 
tion from the homogeneous branch to the inhomogeneous 
branch takes place at a lower value of the energy, and (ii) 
the physical caloric curve presents a region of negative 
specific heat U'(T) < 0. It was noted in |18l26j that these 
anomalies cannot be explained in terms of Lyndon-Bell's 
theory so that other approaches could be considered... We 
here argue that these anomalies can be explained naturally 
by assuming that the QSSs in that region are polytropic 
distributions with an index close to the critical one n c = 1 
(i.e. 7 C = 2 or q c = 3). 

We make the link between the notations of Antoni & 
Ruffo [8] and our notations jQj by setting U — e/4 + 1/2 
and T = 0/2, where U is the energy and T the kinetic 
temperature. This amounts to taking k = 2tt /N and M = 



N in the previous formulae and adding a constant in the 
energy (ground state). For n = 1, we find that the transi- 
tion from the homogeneous states to the inhomogeneous 
states takes place for e c = 1/2, i.e. U c = 5/8 ~ 0.625. 
On the other hand, for indices n* < n < no, the physical 
specific heat U'(T) is negative. In particular, for n = 1, 
we have Chin — ^7 — 11& ~ ~ \- This is at variance with 
Lynden-Bell's prediction. Note, however, that the branch 
of inhomogeneous polytropes exists only above the mini- 
mum energy e = — 7r 2 /16, i.e. above U = — 7r 2 /64 + 1/2 ~ 
0.345 (corresponding to T -> tt 2 /32 ~ 0.308). Below this 
energy, we expect that the system reaches the Lynden-Bell 
distribution. The numerical caloric curve shows a good 
agreement with the Lynden-Bell prediction (equivalent 
to Boltzmann) in that region. However, this agreement 
should be ascertained by explicitly computing the distri- 
bution functions. Similarly, it should be checked whether 
the homogeneous states correspond to Lynden-Bell dis- 
tributions (equivalent to Boltzmann distributions), poly- 
tropic distributions or something else. 

In conclusion, our theory of polytropes qualitatively 
explains the "anomalies" displayed by the caloric curve of 
Antoni & Ruffo [8]. As in previous works |18l26j . we justify 
these polytropic distributions as a result of incomplete vi- 
olent relaxation and lack of ergodicity (inefficient mixing) . 
We stress, however, that polytropes are not universal at- 
tractors in case of non-ergodicity and that other distribu- 
tions could emerge. However, it turns out that polytropes 
play a special role. We have not tried to compare in more 
detail the theory and the observations in order to deter- 
mine the best value of q (it may not be exactly q = 3 
but it happens to be close) . This is left for a future study 
where wc will also determine numerically the distribution 
functions to check whether they can be fitted by q ~ 3 
polytropes. The simulations of Campa et al. |34j go in that 
direction but, since there are differences with the simula- 
tions of Latora et al. [Mj , a detailed study must be per- 
formed. Nevertheless, our approach brings a plausible ex- 
planation of the numerical results of Antoni & Ruffo [5] in 
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terms of Tsallis distributions (polytropcs). Furthermore, 
the apparently selected index q c = 3 is not completely 
arbitrary. First, it leads to distributions with a compact 
support, which is consistent with the phenomenology of 
incomplete mixing. Secondly, it happens to coincide with 
(or be close to) the critical index that appears in the the- 
ory. It is interesting to note that this value was found 
numerically [53] before the theory showed that it has a 
special meaning. All these results have to be confirmed. 



8 Stability analysis 

In the previous sections, we have determined critical 
points of entropy (or free energy) and we have plotted 
the corresponding series of equilibria. For some indices, 
we have found multiple equilibria for the same value of 
E (or K) and, in case of competition, we have compared 
the equilibrium entropy (or the equilibrium free energy) 
of the different solutions in order to select the optimal 
one. However, a more precise analysis should determine 
whether the selected solutions really are entropy maxima 
(or free energy minima) by investigating the sign of the 
second order variations of S (or F). This is a more com- 
plicated problem that will be solved only partially in the 
following sections. 



8.1 Variational principles: general theory 

8.1.1 Microcanonical ensemble 

Basically, we have to solve the maximization problem 

max {S[f] \E[f] = E, M[f] = M} , (162) 

where S[f] is the Tsallis functional pT|) . To solve this 
problem, we can proceed in two steps. First step: we first 
maximize S at fixed E, M and p(r). Since the specifica- 
tion of p(r) determines M and W, this is equivalent to 
maximizing S at fixed Eu n and p{r). Writing the first 
variations as 



SS - p8E kin - / a(r)5 /dv 4 = 0, (163) 



we obtain 



/(r,v) 



mO) 



(q-l)(3v 2 



1/(9-1) 



(164) 



This critical point is the global maximum of en- 
tropy with the previous constraints since S 2 S = 
~ f / f q ~ 2 (3f) 2 drdv < (the constraints are linear in 
/ so their second variations vanish). The function p(r) is 
determined by the density p = J f dv. This leads to the ex- 
pression ([5B]) . From this expression (f5H|) . we can compute 
the pressure p — 4 / fv 2 dv and we obtain the polytropic 
equation of state (|2"7|) . We can now express the energy 
and the entropy as a function of p and p and we obtain 



expressions (j5Tj) and ([5"2]) . These expressions are valid here 
out-of-cquilibrium. Second step: we now have to solve the 
variational problem 



m' AX {S[p]\E[p] 
p 



E\ M[p] = M} , 



(165) 



where S, E and M are given by Eqs. (|5T]). ([52]) and (fT5|) . 
Conclusion: the solution of (|162|) is given by Eq. (|164[) 
[or ([55)1 ] where p is the solution of (|165[) . Therefore, we 
have reduced the study of the initial variational problem 
(|162p for the distribution function /(r,v) to the study of 
a simpler variational problem (|165[) for the density p(r). 
Note that the two problems (|162[) and (|165[) arc equiva- 
lent in the sense that the distribution function /(r, v) is a 
maximum of S[f] at fixed E and M iff the corresponding 
density profile p(r) is a maximum of S[p] at fixed E and 
M. 

The variational problem (|165[) has been studied by 
Chavanis & Sire |46[ in arbitrary dimension of space. The 
first variations return the relation ([55)) . Then, this distri- 
bution is a local maximum of entropy at fixed mass and 
energy iff (see Appendix D of [4*6]): 



1 

- 2 1 



p z 2 



8p6<P dr 



2n 



d(2n — d) j pdr 



<P + -7- I 8 pdr 

2 p 



< 0, (166) 



for all perturbations Sp that do not change mass at first 
order (the conservation of the energy has been taken into 
account in obtaining Eq. (|166[) ). 



8.1.2 Canonical ensemble 

In the canonical ensemble, we have to solve the minimiza- 
tion problem 



min {F[f] I M[f] = M}, 



(167) 



where F = E - TS and T = 1/(3. To solve this problem, 
we can proceed in two steps. First step: we first minimize 
F at fixed M and p(r). Since the specification of p(r) 
determines M, this is equivalent to minimizing F at fixed 
p(r). Writing the first variations as 



(168) 



5F + T a(r)5 / f dv ) dr = 0, 



we obtain 

/(r,v) = 



(q-l)f3v 2 



1/(9-1) 



(169) 



This critical point is the global minimum of free 
energy with the previous constraint since S 2 F = 
2J 1 j f q ~ 2 (5f) 2 drdv > (the constraint is linear in / 
so its second variations vanish). The function p(r) is de- 
termined by the density p = J f dv. This leads to the ex- 
pression (J55J. From this expression ([58]) . we can compute 
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the pressure p — 4 / fv 2 dv and we obtain the polytropic 
equation of state (f2~7|). We can now express the free en- 
ergy as a function of p and p and we obtain expression 
These expressions are valid here out-of-equilibrium. 
Second step: we now have to solve the variational problem 



mm{F[p]\M[p] = M] , 
p 



(170) 



where F and M are given by Eqs. ((55]) and (fT5]) . Conclu- 
sion: the solution of (flUT]) is given by Eq. (fTS9"|) [or {35}] 
where p is the solution of (|170[) . Therefore, we have re- 
duced the study of the initial variational problem (|167[) 
for the distribution function /(r, v) to the study of a sim- 
pler variational problem (|170[) for the density p(r). Note 
that the two problems (|167p and (|170|) are equivalent in 
the sense that the distribution function /(r,v) is a min- 
imum of F[f] at fixed M iff the corresponding density 
profile p(r) is a minimum of F[p] at fixed M . 

The variational problem (|167p has been studied by 
Chavanis & Sire [46j in arbitrary dimensions of space. The 
first variations immediately return the relation (|33[) . Then, 
this distribution is a local minimum of free energy at fixed 
mass iff 



s2p =h 



p- 



dr+l 



5pS<Pdr>0 (171) 



for all perturbations 5p that do not change mass at first 
order. 

Remark 1: from the stability criteria (|166l) and (|171[) . 
we clearly see that canonical stability implies microcanon- 
ical stability (but not the converse). Indeed, since the last 
term in Eq. ()166|) is negative, it is clear that if inequality 
(|171[) is satisfied, then inequality (|166[) is automatically 
satisfied. In general, this is not reciprocal and we may 
have ensembles inequivalence. However, if we consider a 
spatially homogeneous system for which <£> is uniform, the 
last term in Eq. (|166p vanishes (since the mass is con- 
served) and the stability criteria (|166[) and (|171[) coincide. 
Therefore, for spatially homogeneous systems, we have en- 
sembles equivalence. 

Remark 2: the approach developed previously in the 
canonical ensemble can be generalized to any functional 
of the form ^ The minimization problem (|167[) is 
equivalent to (|170[) with the free energy F[p] given by 



F = i / 



P(P') 
P 2 



dp'dr, (172) 



where p{p) is the equation of state associated to C(f). A 
critical point of (|172|) satisfies the condition of hydrostatic 
balance 



Vp = -pV$. 



(173) 



Then, this critical point is a local minimum of free energy 
at fixed mass iff 



6 2 F = 



for all perturbations 5p that do not change mass at first 
order. For the polytropic equation of state (f2"T)) . this re- 
turns the previous results. 



8.2 Application to the HMF model 

Let us now apply the preceding results to the HMF model. 
We start by the canonical ensemble which is simpler in a 
first step. 



8.2.1 Canonical ensemble 

According to the results of Sec. 18.1. 2[ we have to solve 

mm{F\p]\M[p}=M}, (175) 
p 

with 

ttB 2 K 



F[p] = - 



7' 



1 



p 1 dd. 



(176) 



To solve this problem, we can again proceed in two steps. 
We first minimize F[p] at fixed M and B. This gives 



Pi{0) 



1 



-Acos( 



7 



(177) 



where the Lagrange multipliers A and A are determined 
by the constraints M and B through 



.4 



M 



27rJ 7 ,o(A)' 



and 



2ttB _ J 7 ,i(A) 
~kM ~ 7 7 ,o(A) ' 



(178) 



(179) 



This is the global minimum of F with the previous con- 
straints since S 2 F > (this can be deduced from Eq. 
(|17ip by taking S'P = since B is fixed). Then, we can 
express the free energy F as a function of B by writing 
F{B) = F[ Pl ]. This gives 



F(B) = 



ttB 2 



K 



Mi 



x 1 



7-l[27r/ 7i0 (A)]7-i 
7" 



7 



1 2nB\ 



(180) 



where A is related to B by Eq. (fT79")) . Finally, (|T75|t is 
equivalent to 



mm {F(B)}, 



(181) 



in the sense that the solution of (|175p is given by Eqs. 
(fT77|) - (fT79)l where B is the solution of (|T8T|) . For given K 
and M (canonical ensemble), we just have to determine 
the minimum of a function F(B) of the magnetization B. 
The critical points F'(B) = should return the condition 
(|84p of equilibrium. Furthermore, the minima correspond 
to F"(B) > 0. For isothermal distributions, the equations 
take a simple form and we can show [66| through simple 
graphical constructions (by studying the slopes of sim- 
ple curves) that they return the well-known stability re- 
sults. For polytropic distributions, it seems more difficult 
to study (|18ip at a general level. 
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Note that in terms of A, related to the magnetization 
by (|179p . the normalized free energy (|180[) can be written 



/(A) 



A,i(A) 



1-1,0 W 



+ 



1 1 



1 



7- 

7' 



7 



1 r) Ij,o(X) 
£y,o(A) 



7-1 



(182) 



Therefore, for each value of the polytropic temperature 77, 
it suffices to study this function and determine its minima. 
Of course, this can be done easily for each 77 but it seems 
difficult to find a general criterion of stability. It seems 
even difficult to check (in the non-isothermal case) that 
/'(A) = leads to A = x where x is given by Eq. ([54")) . 

However, some analytical results can be obtained for 
the critical index n = 1 (7 = 2). In that case, using Eq. 
(|150[) for A < 2, the expression (|182j) of the free energy 
becomes 



Its derivative is 



(183) 



(184) 



For r\ 7^ rj c = 2, the only critical point /'(A) = is A = 
(homogeneous phase). This is a minimum for rj < r\ c 
(/"(0) > 0) and a maximum for > r\ c (/"(0) < 0). If we 
now consider the critical temperature 77 = rj c = 2, we get 



/(A) = 



1 



(185) 



Therefore, the solutions < A < 2 have the same free 
energy and the inhomogeneous branch is degenerate. 



8.2.2 Microcanonical ensemble 

Let us now consider the microcanonical ensemble. Accord- 
ing to the results of Sec. 18.1.11 we have to solve 



with 



max {S[p]\E[p] = E,M[p} = M}, 
p 



S[ P ] 



n-- )p I pd6, 



E[p] =\Jp 



d6- 



ttB 2 



(186) 



(187) 



(188) 



To solve this problem, we can again proceed in two steps. 
We first maximize S[p] at fixed E, M and B. This gives 
the optimal distribution (|177|) as in the canonical ensem- 
ble. This is the global maximum of S with the previous 
constraints since S 2 S < (this can be deduced from Eq. 
P66|) by taking 8<P = since B is fixed). Then, we can 



express the entropy S and the energy £ as a function of 
B by writing S = S[pi] and E = E[pi\. This gives 



KM 1 



2[27t/ 7 . (A)]7- 



7 



1 2ttB 



ttB 2 



and 



5 



2ri-l 1 



E 



ttB' 



(189) 



(190) 



where A is related to B by Eq. (|179p . The first equation de- 
termines K as a function of B, E and M. Then, S becomes 
a function of B, E and M. Finally, f|l 86[) is equivalent to 



max {S(B)} . 



(191) 



in the sense that the solution of (|186[) is given by Eqs. 
(|T77j) - ([T79l) where B is the solution of (Tl^T|) . For given 
E and M we just have to determine the maximum of a 
function S{B) of the magnetization B. The critical points 
S'{B) = should return the condition (|100[) of equilib- 
rium. Furthermore, the maxima correspond to S"(B) < 0. 
Unfortunately, it seems difficult to study the maximization 
problem (|191[) at a general level, except in the isothermal 
case 7 = 1 [66 . 

Note that in terms of A, related to the magnetization 
by (|179[) , the normalized entropy can be written 



s(A) = - n - - 



1 + 



7-1 \ I-yaW 

7 A Iy,oW 



(192) 



Therefore, for each value of the energy e it suffices to study 
this function and determine its maxima. Of course, this 
can be done easily for each e but it seems difficult to find 
a general criterion of stability. It seems even difficult to 
check (in the non-isothermal case) that s' (A) = leads to 
A = x where x is given by Eq. (|100[) . 

However, some analytical results can be obtained for 
the critical index n = 1 (7 = 2). In that case, using Eq. 
(|150[) for A < 2, the expression (|192[) of the entropy be- 
comes 



s(X) 



1 



We easily obtain 



S '(A) = - 



A 1 



^ 8 



2e- 1 



(193) 



(194) 



The critical points s'(A) = are A = 0, corresponding to 
the homogeneous state, and A 2 = x 2 = 8(1 — 2e), corre- 
sponding to the inhomogeneous solutions (|152p . Then, we 
find that 

s"( X ) = -±=^r y (195) 
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For x < x c — 2, e goes from 1/4 to 1/2. For these values, 
s' (x) < implying that the inhomogeneous solutions A = 
x are entropy maxima (stable). On the other hand, for 
A = 0, we have 



'(0) 



l-2e 
~8e^' 



(196) 



For 1/4 < e < 1/2, we find s"(0) > implying that the 
homogeneous solution A = is an entropy minimum (un- 
stable). 



9 Conclusion 

The nature of quasi stationary states (QSS) in Hamilto- 
nian systems with long-range interactions has created a 
lively debate in the statistical mechanics community |67j . 
Some researchers have argued that these QSSs could not 
be explained in terms of ordinary statistical mechanics and 
that it was necessary to use a new thermodynamics (Tsal- 
lis generalized thermodynamics) [6813] . Other researchers 



have argued, on the contrary, that Tsallis thermodynam- 

This led to a 



ics was unsuccessful to explain QSSs [2412] 
violent polemic. We have since the start adopted an inter- 
mediate position |16l25ll8l26l37j . In principle, the QSSs 
can be explained in terms of Lynden-Bell's theory, i.e. or- 
dinary thermodynamics adapted to the Vlasov equation. 
However, this approach, as usual, assumes ergodicity and 
efficient mixing. There are cases where the system mixes 
well so that Lyndon-Bell's prediction works well [15]. How- 
ever, there are known situations were mixing is not effi- 
cient enough so that other distributions emerge [14134] . 
This is always the case in astrophysics (Lynden-Bell's the- 
ory fails to descri be g alaxies) and in some situations of 2D 
turbulence [28-30]3 Since there is no rigorous theory of 
non-ergodic (partially mixed) systems, we are open to new 
suggestions including the one by Tsallis. 

The idea of Tsallis generalized thermodynamics, as we 
understand it, is to account for non ergodicity and incom- 
plete mixing. It is based on the postulate that, in case of 
incomplete mixing, the system still maximizes an entropy 
but this is not the Lyndcn-Bell entropy. The Lyndcn-Bcll 
entropy is obtained from a combinatorial analysis assum- 
ing that all microstates are equiprobable. This is the main 
postulate of statistical mechanics. Tsallis generalized ther- 
modynamics rejects this equiprobability postulate and ar- 
gues that complex systems "prefer" some regions of phase 
space rather than others. Then, Tsallis form of entropy is 
introduced as a postulate or justified from an axiomatic 
basis. We think that some works remains to be done in 
this area to determine to which situations of incomplete 
mixing Tsallis thermodynamics applies. Tsallis entropy 



24 In the plasma experiment, although Lynden-Bell's predic- 
tion fails to describe the details of the distribution, it provides 
however a fair first order prediction [3D]. The success or failure 
of Lynden-Bell's prediction can be more or less severe depend- 
ing on the initial conditions and it is not possible to draw a 
general conclusion (see the analogy with the meandering course 
of the Mississippi River in [27] ). 



certainly does not describe all types of incomplete mix- 
ing (e.g. galaxies are not stellar polytropes) but it may 
describe some of them. Therefore, Tsallis distributions 
(polytropes) are not universal attractors, but they are not 
useless neither. The results of this paper, and those re- 
ported in |34l37j . suggest that Tsallis distributions can 
be useful to describe QSSs in some situations where the 
Lyndcn-Bcll theory fails. This happens precisely in the 
anomalous region discovered long ago by Antoni & Ruffo 
[8], which stimulated all the activity on the HMF model. 
Therefore, our findings tend to show that the approaches 
of ordinary thermodynamics (Lyndcn-Bell) and general- 
ized thermodynamics (Tsallis) are complementary rather 
than antagonistic (this is often the final outcome of a con- 
flicting situation) [T5] . 

We can also avoid entering in the polemic by tak- 
ing the following simple position [47139125149] . We can 
argue that, in case of incomplete violent relaxation, the 
system reaches a QSS that is a steady solution of the 
Vlasov equation on the coarse-grained scale. This QSS dif- 
fers from the Lynden-Bell prediction but it is not possible 
to predict it as it depends strongly on the dynamics. Of 
course, we must select only nonlinearly dynamically stable 
steady states and this leads to investigate the variational 
problems of the form ((U)-©. These variational problems 
are similar to thermodynamical variational problems^ 
but they lead to dynamical stability criteria, not ther- 
modynamical stability criteria. This dynamical approach 
docs not attempt to explain how the QSSs are selected. 
It just provides formal dynamical stability criteria for a 
large class of distribution functions. In that dynamical ap- 
proach, polytropic distributions just appear as particular 
steady states of the Vlasov equation |47j . By contrast, the 
thermodynamical approach of Lynden-Bell, or its gener- 
alization [15] , seeks to explain how the QSSs are selected 
and predict them from general (yet possibly arguable) con- 
siderations. We think that dynamics and thermodynamics 
are intermingled [47149] and that much work remains to 
be done in this area. Also, the methods of chaos can be 
very useful to explain incomplete relaxation and lack of 
ergodicity [69170] . The subject is certainly not closed. 



A The condition of hydrostatic equilibrium 

For any distribution function of the form / = /(e) where 
e = v 2 /2 + <P(r) is the individual energy, using Eq. (|2"T]) . 
one has 

Vp = iv<P J f'(e)v 2 dv = iv<£ J (j£ ■ v^j dv 



-_V<P / /V v -vdv = -V<2> / fdv = -pV$, 



(197) 



which is the condition of hydrostatic equilibrium ([22 



This is the "Canada-Dry" analogy: in these variational 
problems, S has the color of an entropy, but it is not an en- 
tropy! 
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B Derivation of the polytropic density and C Derivation of the expression (T5TT) of the 
pressure laws entropy 

For n > d/2, the polytropic DF can we written In the expression (fTTj) of the entropy, we have to evaluate 

J f q dv. Let us first consider the case q > 1. Using Eq. 
/ = A(e m - t) n ~ d/2 . (198) (J3HD and defining x = v 2 /[2(n + l)Kp^ n ], we get 

The density and the pressure can be expressed as f ra i 1 ^ n n ± z 2 r/ 1 N r ^ d / 2 



P = AS d Q ($), P = -AS d Q 2 {$), (I99) ,i 



/ ' f q dv = ^p r S d ? i r [(n + l)K}° 



with 



{l-x) n+1 -ix^ dx. (208) 



The integral can be evaluated in terms of r functions and 
' v ' /s+d-i t , /onn'i we obtain 



Q fc =/ ^e m -<?-y) w^*- 1 ^. (200) 

Setting a? = w 2 /[2(e m — #)], we obtain 



| f<dv = ±prS d 2*? [(n + l)Kf 2 



r(d/2)r(2-d/2 + n) , . 

Qfe = 2 (^-2)/ 2(em _ #) n + */a V( 2 + n) ' (2 ° 9) 

x J (1 — x ) ^ x^ ~ + ^ dx. (201) Using the property .T(x + 1) = xT(x) and the expression 

(139)) of Z, the foregoing equation reduces to 



The integral can be expressed in terms of Gamma func- 
tions leading to f p dv = ^ p i 1 ~ d / 2 + n _ (210) 



r{{d + k)/2)r{\- d/2 + n) 



(202) 



Zl- 1 ^ 1 + 71 

Using the equation of state (j27|) . we obtain 



r<1 + ' /2 + n) /^-T^fe 

Then, the density and the pressure can be expressed in J V ' ; 

terms of the potential <P as in Eqs. (E3D - dMJ) ■ _ Combining Eqs. @ and ® , it is easy to establish that 
lor n < — 1, the polytropic Dl can we written 

f = A(e + e) n - d / 2 . (203) (n + \)KZ^ = (212) 

The density and the pressure can be expressed as Substituting this identity in Eq. (J2TT]), we find that 



P 

with 



AS d R (<P), p = - d AS d R 2 {$), (204) J f 9 dv= { 1 _* + A «-i p . 



(213) 



+oo / 2\ "—d/2 



Recalling the definition of A after Eq. (|T9")) , we obtain 
i?fe = / (£0+^+^-) u fc+d_x du. (205) f ( d \ o-l 

Setting a; = w 2 /[2(eo + #)], we obtain 



q 

After simplification, this yields the simple result 



R k = 2 {k+d - 2)/2 (e + <P) n+k / 2 

+00 / f q dv = /3p. (215) 

x / (l + x) n - d ' 2 x^ +d - 2 ^ 2 dx. (206) 



Substituting this identity in Eq. ([TTj) . we obtain Eq. (|5T1) . 
The integral can be expressed in terms of Gamma func- The case q < 1 can be treated similarly and yields the 
tions leading to same final result. 

Rk = 2 (*+<*-2)/2( eo+( p)»+fc/2 
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